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abstract 


1 aN  Identify  N Meek 


is  report  presents  the  development  and  exposition  of  a non-Gaussian  aerosol 
dispersion  model  and,  the  effect  of  the  aerosol  cloud  on  electrooptlcal  systems 
Explicitly,  algorithms  are  provided  to  calculate  both  the  transmitted  and 
backscattered  (glare)  radiation  components. 

Finally,  a validation  of  the  model  was  performed  with  good  to  excellent  agree— 

• This  validation  Included  the  space— time  history  of  cloud  geometrv  and 
aerosol  concentration. 
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SUMMARY 


An  approach  to  the  effect  of  smoke  or  obscuring  aerosols,  on  battle- 
field electrooptlcal  systems  Is  presented  here;  that  Is,  a first  gene- 
ration model  Is  constructed  to  be  consistent  with  an  earlier  developed 
GRC  target  acquisition  model  for  CARMONETTE  or  any  combat  simulation 
model.  The  attack  taken  here  is  to  consider  smoke  as  a transient  modi- 
fication to  the  natural  terrain  height,  l.e.,  the  spatial  coordinate 
(x,y,z)  become  [x,y,z+Z(t)]. 

The  approach  philosophy  was  to  Initially  generate  a simple  semi- 
emplrical  model  that  Includes  all  the  necessary  ingredients  to  permit, 
if  necessary,  the  generation  of  a model  of  greater  sophistication  and 
accuracy.  The  model  is  not  the  classical  approach  (Gaussian);  however, 
the  model  validates  extremely  well  with  field  tests. 

In  Volume  1 a technical  discussion  is  presented  and  includes  smoke 
generation  mechanisms,  smoke  dispersal,  radiation  propagation  and  vali- 
dation tests  of  the  proposed  smoke  model.  In  Volume  2 we  present  the 
model  structure  in  a step-by-step  procedure.  / 
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BACKGROUND 


This  report  is  a formalization  and  extension  of  an  earlier  working 
paper1  prepared  for  the  U.S.  Army  Night  Vision  Laboratory  on  a model 
to  describe  the  smoke  cloud,  aerosol  concentration  and  optical  proper- 
ties of  the  cloud  where  the  smoke  is  generated  by  various  smoke  munitions. 

The  diffusion  of  smoke  and  other  pollutants  have  been  studied 
for  the  past  forty  years,  both  theoretically  and  experimentally.  Since 
the  primary  transport  mechanism  is  atmospheric  turbulence,  the  overall  activity 
has  been  directed  in  solving  the  turbulent  diffusion  equation  and  measuring 
the  cloud  parameters  in  terms  of  solution(s)  to  the  latter  equation,  i.e, 
in  terms  of  mean  values  and  their  variances.  The  difficulty  with  the 
application  of  this  approach  is  that  (1)  the  use  of  the  "turbulent" 
diffusion  equation  to  find  the  mean  concentration  and  mean  cloud  dimensions 
is  normally  based  on  3 minute  averaging  periods  or  longer  and  (2)  the 
boundary  condition  for  the  near  ground  source  is  inadequate.  Of  the  two 
issues  the  former  may  be  more  critical  since  in  most  cases  of  interest 
the  smoke  generation  is  over  in  at  most  several  minutes  and  fluctuation 
statistics  are  not  satisfied  by  the  suggested  sampling  time. 

Also,  there  is  some  question,  from  purely  a mathematical  standpoint 
whether  the  Sutton-Calder  solution  (parabolic  diffusion  equation)  is 
appropriate  In  place  of  a hyperbolic  one  (Monin-Oboukov). 

Hence,  for  the  present  purpose,  namely  to  generate  a simple  model 
that  can  be  utilized  in  a combat  simulation  model  an  empirical  approach 


u 
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is  adopted  throughout  its  formulation  so  that  the  degradation  effects 
of  smoke  on  military  systems  can  be  readily  calculated.  -• 

This  volume  Is  basically  divided  Into  two  parts:  Part  1 will  be 
devoted  to  the  concepts  related  to  the  model;  and  Part  2 will  deal  with 
the  model  proper.  Its  formulation  and  several  validation  tests.  Whenever 
ancillary  data  are  required  these  will  be  provided  In  one  of  several 
appendices. 

The  writer  wishes  to  express  his  appreciation  to  A.  Linz  and  R. 

Bergeman,  Night  Vision  Laboratory, for  their  encouragement  during  the 
initial  study  phase,  Seymour  Berler,  DIA,  for  his  interest  and  agency 
support,  and  L.  Obert,  NVL  COTR,  for  review  of  this  report. 
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1 . INTRODUCTION 


1.1  OBJECTIVE 

The  primary  goal  is  to  develop  a methodology  for  the  introduction  of 
a smoke  model  into  CARMONETTE  that  will  handle  both  passive  battlefield 
and  active  battlefield  electrooptical  systems,  visual  and  infrared. 

This  includes  target  acquisition  devices,  missile  seekers  and  laser 
systems. 

1.2  APPROACH 

We  consider  here  smoke  as  an  artificially  controlled,  obscuring 
aerosol  cloud  which  can  reduce  the  transmittance  along  the  line  of  sight 
between  the  observer  and  the  target.  In  addition,  backscatter  or  self 
emission  from  the  cloud  may  reduce  the  apparent  contrast  at  the  sensor 
input. 

Thus  for  passive  and  active  systems  operating  in  visible  (photopic), 
near  infrared  (NIR),  mid  infrared,  or  far  infrared  (LWIR)  wavelengths 
we  must  consider  the  wavelength  dependence  of  the  optical  properties  of 
smokes.  It  should  be  noted  that  there  are  a variety  of  smokes*  with 
different  optical  characteristics.  In  order  to  establish  the  macroscopic 
effects  of  a smoke  cloud,  the  spatial  distribution  of  the  aerosols  is 
required  as  a function  of  time.  The  latter  depends  in  a complex  way  on 
the  terrain,  meteorological  conditions,  character  of  smoke  munitions 
employed,  etc.  The  subject  Is  treated  in  a statistical  fashion. 

Prom  a combat  simulation  viewpoint,  smoke  in  the  line  of  sight  is 
equivalent  to  a time  dependent  loss  of  intervisibility.  The  approach 

*Smokes  and  obscuring  aerosols  are  used  in  an  equivalent  sense. 


suggested  here  is  to  eAumlne  the  proposition  that  a smoke  cloud  or  screen 
be  treated  as  a time  dependent  terrain  feature,  that  is  the  incremental 
Increase  in  the  spatial  coordinates,  see  schematic  diagram. 


The  terrain  coordinates  at  "A"  are  (x,y,z)  and  with  smoke  became  (x,y, 
z(t) ) ; i.e.,  not  unlike  vegetation. 

The  advantages  of  this  approach  are  (1)  inputting  a time  dependence 
for  smoke  effects  Is  compatible  with  the  CARMONETTE  war  game,  (2)  the 
Integration  of  smoke  with  terrain  may  be  a significant  aid  for  "intel- 
ligence preparation  for  the  battlefield"  and  (3)  the  concept  handles 
intervlslbllity  for  various  weapons,  i.e.,  ground-ground,  ground-air, 
air-ground,  air-air,  and  target  acquisition  systems. 

Briefly  then,  we  require  the  formulation  of  the  following: 

(1)  The  general  properties  of  smoke  munitions 

(2)  Temporal -spatial  history  of  a smoke  plume  an  i 

(3)  The  radiation  transfer  properties  of  the  plume. 
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The  amount  of  data  on  each  of  the  above  is  limited  particularly 
as  the  degree  of  detail  increases:  (1)  for  hygroscopic  smokes,  the 
quantity  of  energy  released  is  required  to  determine  the  smoke  plume 
temperature  history;  (2)  a valid  description  of  the  smoke  plume  is  vital 
for  the  prediction  of  thermal  systems  and  the  turbulence  effects  on 
laser  systems;  and  (3)  the  particulate  distribution  will  affect  the 
radiation  transfer  properties. 

In  the  latter  instance,  the  transmittance  and  glare  may  be  esti- 
mated; however,  the  inhomogeneities  generated  by  random  distributed 
smoke  munitions  can  seriously  affect  thermal  systems  operations  beyond 
pure  transmittance  and,  must  be  taken  into  account. 


2.  BASIC  CONCEPTS 


2.0  INTRODUCTION 

We  will  discuss  here  a number  of  basic  concepts  needed  to  sup- 
port the  smoke  model  that  is  to  be  presented  later  and,  of  course,  its 
limitation.  These  concepts  cover  sources,  cloud  dynamics,  and  optical 
properties. 

2.1  SOURCES 

The  smoke  munitions  used  by  various  countries  can  be  categorized 
as  (1)  non-exothermic  aerosol  generators  and  (2)  exothermic  aerosol 
generators.  A primary  example  of  the  latter  is  white  phosphorus.  Next 
we  can  divide  these  into  instantaneous  (short  burn  time)  and  finite  burn 
time.*  Finally,  we  may  have  "point"  sources  and  finite  line  sources. 

The  latter  is  accomplished  by  detonating  several  munitions  or  submunitions 
in  a linear  array,  a linear  array  of  smoke  pots  or,  a linear  line  of 
spray.  For  completeness  several  standard  munitions  are  given  in 
Appendix  A. 

2.2  YIELD 

There  are  a variety  of  chemical  compositions  that  may  be  used  for 
smoke  generation  namely  white  and  red  phosphorus  (WP/RP),  hexachlorethane 
+ Zn  (HC),  fog/dlesel  oil,  TiCl^(FM) , chlorosulfonic  acid  + S03(FS), 
sulfuric  acid  + SO^Oleum),  and  anthracene  mixtures.  The  products  may 

*We  wish  to  make  a distinction  between  reaction  time  to  convert  the 
munition  fill  into  smoke  and  smoke  cloud  persistence.  Here  we  utilize 
the  time  to  convert  primary  material  into  smoke. 
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range  *rom  mists  of  oil  droplets,  metallic  chlorides  (ZnCl2,  TiCl^ 
NH4CI,  etc.)  to  acid  droplets  like  sulfuric  acid  or  phosphoric  acid 
(H3PO4).  The  interaction  of  these  aerosol  products  with  atmospheric 
moisture  varies  considerably;  that  is,  they  may  absorb  moisture  and  grow 
or  react  chemically  like  phosphorus  pentoxide  to  form  a phosphoric  acid 
and,  thereby  increase  in  weight.  This  property  has  been  parameterized 
as  the  smoke  yield  where  the  yield  is  defined  by  the  expression 


v _ weight  of  smoke  generated  n* 
~ weight  of  smoke  fill 


(2.1) 


weight 

It  should  be  noted  that  the  significance  lies  in  the  fact  smoke  concen- 
tration is  given  by  gms  of  smoke/unit  volume. 

Johnson  and  Forney  (see  Appendix  B)  show  that  the  yield  is 


dependent  on  the  relative  humidity  (RH)  for  those  aerosols 
affected  by  atmospheric  moisture.  The  expressions  for  several  smokes 
are  given  below  where  the  constant  applies  for  RH  £ bracket  value: 

Y(fog  oil)  =1 

Y(WP)  = .003  (RH  - 40)1,67  + 2.9 

85  (2.2) 

Y(ZnCl2)  = .051  (RH  - 50)‘a3  + 1 

Y(A1C13/H2S04)  = .016  (RH  - 20)1'25  + 1.4 

The  use  of  Y - 1 for  fog  oil  Is  appropriate  with  the  possible  exception 
at  high  relative  humidity  since  the  oil  may  adsorbed  moisture. 

Clearly  the  yield  Is  related  to  the  hygroscopic  properties  of  the  con- 
densation nuclei  produced  by  the  reactions  of  the  smoke  material.  The 
nuclei  have  radii  between  0.2  - 0.5  pm  with  a mode  of  ~0.25  pm  and, 
subsequently  the  smoke  particles  will  grow  to  a mean  value  of  -1.2  pm. 
*n  is  munition  efficiency  and  is  assumed  to  be  equal  to  1. 
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The  use  of  relative  humidity  rather  than  absolute  humidity  in  eq. 

o 

(2.2)  is  further  substituted  in  measurements  by  Novikov  et  al.  In  Fig.  1 
we  present  the  ratio  of  liquid  mass/solid  mass  vs  relative  humidity  and 
absolute  water  content  for  Fe4Cl/AlCl 3/phosphoric  acid  aerosols.  From  Fig.  1 
we  can  deduce  that  the  yield  is  better  represented  as  a function  of  relative 
humidity  since  the  ordinate  is  proportional  the  yield  minus  one.  The  rela- 
tive humidity  effect  should  be  reflected  in  the  optical  properties  of  the 
smoke  cloud. 

2.3  SMOKE  CLOUD 

To  familiarize  the  reader  we  will  present  here  an  abbreviated 
overview  of  concepts  related  to  the  diffusion  of  smoke  and  other 
pollutants  by  the  atmosphere.  For  more  general  discussion  we  suggest 
refs  3,  4,  5. 

2.3.1  Stability  Parameter 

The  behavior  of  the  smoke  cloud  is  determined  primarily  by  the 
turbulent  properties  of  the  atmosphere  within  which  the  smoke  is  being 
generated  and  therefore,  is  with  the  exception  of  phosphorus  a momentum 
exchange  process.  In  the  case  of  phosphorus  or  other  strong  exothermic 
smokes,  the  rate  of  heat  generated  by  the  associated  exothermic  reactions 
is  the  initial  dominant  mechanism,  i.e.,  for  about  200  secs  or  less. 

For  the  present  we  will  discuss  the  momentum  dominant  smoke  cloud. 

2 

Novikov,  Yu.  I.  et  al.  Advances  in  Aerosol  Physic:-  No.  2,  90,  Israel 
Program  for  Scientific  Translations,  1970. 

3 

Pasquill,  F.,  Atmospheric  Diffusion,  Van  Nostrand,  London  1962. 

4Briggs,  G.  A.,  Plume  Rise,  ESSA,  Oak  Ridge,  Tenn,  1969  (TID-25075). 
Meteorology  and  Atomic  Energy  (D.  H.  Slade,  ed),  AEC,  1968  (TID-24190). 
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Fig.  1 (Liquid/Sol Id)  Aerosol  Mass  vs  Relatlv 
Humidity  (S),  Absolute  Uater  Density 


a) 


w 


. 1 75 ( RH ) - n.62 

.18  (RH)  - 12.72 
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To  categorize  the  atmospheric  capability  to  diffuse  a smoke  cloud 
Pasquill,  British  Meteorological  Office,  introduced  a stability  classi- 
fication which  depends  on  net  radiation  flux  near  the  surface  and  wind 
speed.  Clearly  the  former  is  dependent  upon  diurnal  variations,  cloudi- 
ness, soil  moisture,  etc.  The  three  basic  stability  classes  are  (1) 
unstable,  (2)  neutral  and  (3)  stable  which  can  be  further  refined  to 
several  subclasses  of  unstable  and  stable  conditions.  Neutral  conditions 
exist  with  cloudy  skies  or  high  wind  speeds.  Hence  the  Pasquill  classi- 
fication system  runs  from  A (extremely  unstable)  to  G (extremely  stable) 
with  D * neutral,  see  Table  1 and  Appendix  C. 

Table  1 

Pasquill  Stability  Categories 


Night 


Surface  wind 
speed  at  10  m 
(m  sec*') 

Insolation 

Thinly  overcast 
or  »4/8 
low  cloud 

*3/8  cloud 

Strong; 

Moderate 

Slight 

2 

A 

A-B 

B 

-J 

2-3 

A-B 

B 

C 

E 

F 

3-5 

B 

B-C 

C 

0 

E 

5-6 

C 

C-D 

D 

0 

D 

>6 

c 

0 

0 

D 

D 

It  should  be  noted  that  the  classical  stability  definition  is  given  by 
by  the  atmospheric  lapse  rate  relative  to  the  adiabatic  lapse  rate,  l.e.. 

Potential  temperature,  0(T)  * AT  + r 
where  AT  = T(10m)  - T(0.5m)  and  r * .098cC/10m.  When  9(T)  * 0 the 
condition  is  stable,  neutral  and  unstable  respectively.  An  atmospheric 


stability  parameter  is  given  by 

g . a ii 

5 T 5z 

where  g and  T are  gravitational  acceleration  (m/sec2)  and  ambient  ground 
temperature  (°K)  respectively.  This  parameter  Is  used  later  in  eq.  (3.15d). 

From  ref  (6)  we  reproduce  for  completeness  Table  2 wherein  a contin- 
uous representation  of  the  Pasqulll  classification  is  given.  Here  Class 
A * -3,  Class  G * 3 and  the  neutral  condition  » 0. 


Table  2 

Key  Used  to  Estimate  the  Continuous  Stability  Classification* 


Transient  period 


Night 


Mean 

Kind 

Cm  «rc"*i 

Incoming  solar  radiation 

Strung  Moderate  Slight 

-3.3 3.0 

-3.0 2.2 

-2.3 2.0 

z-.l 

-JO- — 2.2 

-2.2- -2.0 

—2.0 1.0 

1-5 

-2.2 1.3 

-2.0 1.0 

-1.0— 0.3 

K<i 

-1.3 1.0 

— 1.0— OJ 

-0.3 OJ 

— 1.0 0-3 

— 0J-— 0.1 

— OJ—  0 

-0J 

-0.1 

0 

Tbkiir 
overcast 
or  4/8 

Day—  -Night  low  clnuil  Cloud  (3/8 


-U— 0.5 

oj~ij 

1J 

2J 

-1.0 OJ 

QJ~1.0 

1J— 0.6 

2.3— 1.6 

—1.0 — OJ 

0J~0.5 

ao^oj 

1.6— OJ 

-0.4 0.2 

0.2-0.4 

OJ'-O.l 

0.3-.OJ 

— OJ— 0.1 

0.1^42 

0.1~0 

0J-0.1 

• 0 

0 

0.1 

* : I)  Ser  Turner  (1949)  for  overall  ^tplanatioos. 


To  appreciate  the  statistical  variations  in  Pasqulll  Categories 
we  provide  in  Table  3,  the  values  for  291  stations  throughout  the  U.S. 
and  several  stations  In  Hawaii  and  Puerto  Rico  where  the  column  legends 
are  self  explanatory. 

^Shlr,  C.C.  and  Shieh,  L.  F.,  Jour.  Appl.  Met.,  Mar  1974,  p.  189. 
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Table  3 


Statistical  Summary  of  Percent  Frequencies 
of  Pasquill  Stability  Categories 


• 

A 

8 

C 

Unacabla 

(A+a+C) 

Naucral 

(D) 

Stable 

(E+F+C) 

Annual :Me an 

1.0 

6.6 

11.2 

18.8 

47.7 

33.4 

Mod* 

1 

4 

11 

19 

34 

28 

Scd.  Dav. 

1.1 

3.0 

2.2 

3.2 

12.2 

7.3 

Ranga 

0-3 

1-16 

3-19 

6-33 

13-77 

13-54 

314  acna. 

Wlncar:Maan 

0.0 

2.6 

7.1 

9.7 

57.5 

32.7 

Hoda 

0 

1 

6 

7 

66 

26,30,32 

Scd.  Dav. 

0.2 

2.3 

2.9 

3.1 

15. i) 

10.4 

Ranga 

0-1 

0-14 

2-18 

2-28 

14-88 

10-66 

291  atss. 

Spring :Ma  an 

1.0 

6.2 

11.2 

18.4 

32.2 

29.3 

Hoda 

1 

4,6 

11 

14 

62 

23 

Scd.  Dav. 

1.2 

3.0 

2.3 

3.7 

12.6 

7.3 

Ranga 

0-6 

1-16 

3-21 

6-36 

13-79 

12-33 

291  acna. 

SuaaartHaan 

2.6 

11.3 

16.4 

30.2 

36.1 

33.7 

Hoda 

1 

9 

16 

34 

32 

38 

SCd.  Dav. 

2.3 

4.0 

2.3 

3.8 

12.5 

7.7 

Ranga 

0-13 

2-21 

9-29 

13-49 

12-73 

11-49 

291  acna. 

Aucian:Maan 

0.3 

6.2 

10.2. 

16.9 

44.3 

38.3 

Hoda 

0 

6 

9 

18 

38 

42 

Scd.  Dav. 

0.8 

3.7 

2.6 

6.2 

13.9 

8.1 

Ranga 

0-3 

1-19 

3-19 

4-34 

9-73 

18-57 

291  acna. 

Table  4 

Prevalence  of  Pasquill  Stability  Categories 
over  the  Netherlands* 


Ooy-WtfMr 


Doy-Somq 




Nigm 

cftt  Motors 


V 'Wind  toMO- '"on 
N • Tara  cou<jm«*»  * 
* > Cdscjrtd 


hs 


3ay-?«i 


♦Reproduced  from  G.  H.  Strom,  "Atmospheric  Dispersion  of  Stack  Effluents" 
p.  257;  Air  Pollution  (A.C.  Stern,  Ed)  2nd  Ed.  Academic  Press,  1968. 
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The  above  discrete  or  continuous  Pasqulll  system  suffers  from  the 
absence  of  surface  roughness  dependence.  Golder7  established  a relation 
between  the  Monln-Obukov  length,  L,  and  Pasqulll  stability  namely, 

1/L  -±  [d  In  (1.2  + 10/Zo)]2  10  f(s),  (2.3) 

f(s)  ■ -a/(l  + b|s|c) 

s » stability  class  defined  in  Table  2. 
a = 4,  b * 1.3,  c = 0.85,  d = 0.21659  and 
Z0  * is  the  surface  roughness  in  meters. 

The  surface  roughness,  Z0,  varies  from  ~10“3m  (sand)  to  -.10  m (tall 
vegetation).  Kung  and  Lettau  have  found  an  empirical  relation  between 
Z0  and  plant  height,  hc,  l.e. , 

log  Z0  * -1.24  + 1.19  log  hc  (2.4) 

Several  values  of  Z0,  hc  are  given  in  Table  5. 


Table  5.  Profile  Parameters  of  Some  Crops  and  Buildup  Areas 


Crop 

Z0(m) 

hc(m) 

Plane,  snow  covered 

.0049 

.03 

Grassy  surface 

.0173 

.10 

Low  grass 

.032 

.20 

High  grass 

.039 

.30 

Wheat 

.04 

1.3 

Suburban  area 

.4 

Urban 

6 

24 

1 

li  1 

Golder,  0.  "Relations  Among  Stability  Parameters 

In  the  Surface  Layer," 

Boundary  Layer  Meteorology  3 (1972),  p 47. 

o 
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— 

and 

where 


The  Monin-Obukov  length  is  defined  by  the  relationship 


where 


friction  velocity  («tl)  (m/sec) 
acceleration  of  gravity  (m/sec2) 
mean  temperature  (z  s 10m) 
mean  wind  speed  (z s 10m) 


The  effect  of  surface  roughness  is  to  reduce  in  some  instances  the 
stability  class;  i.e.,  say  B->C  or  D-*£,  see  Fig.  2. 


Fig.  2.  1/t  as  a Function  of  Pasquill  Classes  and  Z0 


Evidently  a factor  that  plays  a major  role  in  stability  is  the 
standard  deviation  of  vertical  component  of  turbulence,  ow.  This 
can  be  represented  by  an  approximate  expression, 

crw  (cm/sec)  s 12.3  (u1  - 1) 

where  u^  * wind  speed  (m/sec)  at  1 meter  altitude. 

The  inclination  of  the  wind  speed  is  defined  by  (o^IT)  and, 
therefore, 

fw  a 1 2. 3(u-|  - 1) 

IT  5= • 

For  example,  if  aw  * .5  m/sec  and  IT  = 5 m/sec  then  the  ratio  =0.1 
radians  which  corresponds  to  a slightly  stab1e-*neutral  stability 
condition.  That  is,  for  ow/TT  >0.25  radians  and  IT  s 4 m/sec  an 
unstable  condition  exists,  see  p.  66,  ref  3. 

Gustiness  is  defined  by  where  w'  is  the  fluctuations  in 

the  vertical  velocity  and, is  ~ a^TT  . The  vertical  perturbations  can 
be  estimated  therefrom  provided  u^(t),  the  time  dependent  horizontal 
velocity  at  one  (1)  meter  is  known.  This  may  be  applicable  to  a 
non-exothermic  smoke. 

For  the  present  work  we  adopt  the  simple  stability  definitions 
given  In  Table  1. 

2.3.2  Solution  to  Diffusion  Equation 

We  now  present  a brief  discussion  of  the  Gaussian  solutions  to  the 

diffusion  equation  so  that  these  may  be  compared  to  the  anpirlcal  diffusion 
expressions  used  here. 


16 


The  notation  for  the  following  discussion  is  given  below  as  follows: 

x = direction  of  mean  wind  velocity  vector 
y = perpendicular  to  x 
z = perpendicular  to  plane  (x-y) 

¥ = time  mean  wind  speed  in  x (m/sec) 

V = time  mean  wind  speed  in  y (m/sec) 

W = time  mean  wind  speed  in  z (m/sec) 

t = time  (sec) 

V = volume  (m3) 

Qs  = smoke  munition  fill  weight  (gm) 

Q * smoke  weight  (gm) 

C * concentration  (g/m3) 

X * plume/cloud  length  (m) 

V = plume/cloud  width  (m) 

Z * plume/cloud  height  (m) 

Ki  * turbulent  coefficient  (cm2/sec)  (i  3 x,y,z) 

(a)  Instantaneous  Sources* 

Consider  a smoke  source  at  the  origin  of  a rectangular 
coordinate  system  (0,0,0)  at  time,  t * 0,  then  with  a wind  in  the  x-direction 
having  a wind  mean  speed,  ¥,  the  mean  concentration  of  smoke  will  be 
defined  by 

C(x,y,z,t)  = 8?iFKt)'3/2  exp  -( [x  - ut] 2 + y2  + z2/4Kt)  (2.6) 

*See  discussion  in  Section  2.5. 
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Similarly  for  an  instantaneous  infinite  line  source  we  can  write 
C(x,z,t)  * 4?Kt"  exp  ’ Kx-ITt)2  + z2/4Kt]  (2.7) 

Clearly,  these  solutions  are  travelling  Gaussian  distributions  with 
increasing  entrainment  of  air  to  dilute  the  smoke  concentration.  Further, 
for  the  infinite  line  source  the  diffusion  is  Independent  of  "y". 
Graphically,  the  result  is  depicted  in  Fig.  3. 


Fig.  3.  Schematic  of  Gaussian  Solutions 


(b)  Continuous  Sources 

For  the  point  source  one  obtains  an  expression  of  the  time-mean 
concentration  in  the  form 


C » 


_JL 


4irKr 


[1  -*( 


) 1 


(2.8) 
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where  Q is  the  release  of  smoke  in  (gm/sec)  and  <f>  is  the  normal 
probability  integral  function.  When  a wind  is  introduced  then 


In  the  above  solutions,  turbulence  was  considered  to  be  isotropic,  i.e., 

K = KX  = Ky  * K2,  if  we  assume  no  downward  transport  at  surface,  Kzdr/dz-»-0 
and  conservation  of  downwind  flux.we  obtain  for  a continuous  point  source. 

Z-  (2-10) 

For  a continuous  line  source  of  infinite  length  with  constant  mean  wind, 
tr,  perpendicular  to  y-axis  and  constant  Kz,  the  concentration  is  found 
to  be 

r<x-z)  = exP  -t  ! (2-11) 

where  Kz  can  be  found  from 

u * u 1 (z/z-j  )p*  and 
Kz  = K, 

Subscript  1 refers  to  the  value  at  1 meter.  Further  the  exponent 
"p"  has  values  0<p<l ; for  strong  lapse  p = 0.1,  p~0.8  for  extreme 
stability  and  p = 1/7  for  neutral  stability  and  smooth  ground  surface. 

* This  is  called  the  power  profile  law. 

**This  is  called  the  conjugate  profile  law. 
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(2.12) 


Calder  proposed  a definition  of  u in  a alternate  form  namely, 
u - u*  q (Z/Z0)“ 

where  Z0,  q,  o are  constants  and  u*  is  the  friction  velocity.  Here 
the  value  of  "Z"  is  the  mean  cloud  height  which  is  to  be  established  from 
preliminary  observational  estimates  and  is  required  before  q and  a may  be 
evaluated. 

The  eddy  diffusion  coefficient  Is  defined  by  the  expression 

Kz  = k u*  Zq  (Z/Z0)6  (2.13) 

where  k is  the  Karman  constant  = 0.4  and  Z0  is  the  roughness  length  in 
centimeters. 

Numerical  solutions  to  the  diffusion  equation  are  presented  in 
graphical  form  In  ref  8.  In  addition,  there  is  a thorough  discussion  of 
the  classical  Sutton  treatment  which  is  based  upon  fluctuation  theory  and 
his  solutions  reduce  to  the  above  if  K/TT  = constant. 

2.3.3  Observations  on  Diffusion  Solution 

Irrespective  of  the  origin  of  the  solutions,  i.e.,  (1)  Roberts  - 
presented  above,  (2)  Sutton,  Calder,  etc.,  they  all  suffer  from  major 
deficiencies,  (a)  they  are  steady  state  solutions;  and  (b)  the  cloud- 
ground  plane  interaction  is  inadequate.  For  example.  Shir  (6)  points  out 
that  the  effects  of  surface  roughness  can  Influence  the  concentration 
distribution  significantly.  Similarly,  Hildebrand  (9)  states  as  recently 
as  1977  "there  has  not  yet  been  a systematic  verification  of  the  Gaussian 

— 

Haltlner,  G.  J.  and  Martin,  F.L.,  Dynamical  and  Physical  Meteorology, 
McGraw-Hill,  New  York,  1957. 

^Hildebrand,  G.,  J.A.M.  1977. 
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solution  of  the  equation  for  differing  atmospheric  conditions."  Never- 
theless, for  sustained  operations,  i.e.,  where  smoke  pots  or  generators 
are  utilized  whose  duration  > 30  mins,  use  of  the  Calder-Sutton  equation 
may  be  a fair  approximation. 

2.3.4  Thermal  Effects 

In  the  above  discussion  we  have  not  indicated  any  dependence  on 
thermal  effects  due  to  munition  material.  The  magnitude  of  the  effects 
will  vary  from  minor  (smoke  pots)  to  major  (white  phosphorus)  and  thereby 
alter,  to  varying  degrees,  the  smoke  cloud  dynamics.  We  will  now  examine 
several  effects. 

The  temperature  differential,  Tsmo|<e"Tam|:^ent  al.r>determines  the 
buoyancy  or  rise  and,  therefore  the  rise  will  be  less  for  an  inversion 
(stable)  than  for  lapse  (unstable)  condition. 

For  non-explosive  or,  weak  to  non-exothermic  reactive  smoke  muni- 
tions (material)  the  resultant  initial  rise  is  insignificant  about  1-2 
meters  subsequently  the  cloud  dynamics  is  purely  mechanical,  i.e.,  momentum 
exchange  only. 

For  explosive  munitions  where  the  HE^  will  release  about  5x10^ 
watts/#  the  rise  for  a standard  shell  (105/155nm)  may  be  < 5 meters.  Again 
this  may  be  considered  as  an  initial  rise  upon  which  is  superimposed  the 
normal  rise.  When  a white  phosphorus  round  is  utilized  the  strong  exother- 
mic reaction  will  sustain  the  rise  since  the  heat  content  stays  for  a long 
period  and,  therefore,  strongly  interacts  with  the  mechanical  forces. 

This  topic  has  received  modest  attention  with  respect  to  modelling 
of  smoke  generated  by  explosive  munitions. 


Lastly,  the  ABC-M3A3  pulse  jet  smoke  generator  is  a complex  device  because 
it  emits  about  sixty  (60)  puffs  of  hot  oil  per  second.  To  a first  approxi- 
mation one  may  assert  that  a continuous  source  representation  is  adequate 
and  valid.  The  excess  exhaust  temperature  is  expended  by  the  expanding 
jet  and,  therefore,  the  cloud  rise  is  not  unlike  any  smoke  pot. 

2.4  OPTICAL  PROPERTIES 

The  optical  properties  of  the  different  product  aerosols  and  obscur- 
ants vary  with  wavelength  and  composition,  see  Appendix  E.  As  noted  earlier, 
the  metal  chlorides  and  phosphorus  compounds  are  hygroscopic  and,  there- 
fore, behave  differently*  In  the  visible,  near  Infrared  and  Infrared.  Lab- 
oratory or  field  observations  performed  at  low  relative  humidity  < 40%  may 
be  entirely  Inaccurate  if  corrections  are  not  made  when  the  relative  humidity 
exceeds  50%.  The  errors  arise  from  two  sources  (1)  the  particle  size  dis- 
tribution will  change;  and  (2)  the  optical  properties  will  change  similar 
to  those  observed  for  ambient  atmospheric  aerosols.  The  effect  of  the 
distribution  function  Is  primarily  of  Interest  to  a detail  radiation  trans- 
fer calculation  whereas  the  other  error  can  lead  to  underestimates  of 
transmittance  or  overestimates  of  backscatter. 

2.5  OBSERVATIONS 

From  the  above  discussion  the  formulation  of  an  analytical  model 
for  smoke  clouds  generated  by  munitions  and  the  ability  to  establish  the 
degradation  Induced  In  military  optical  systems  Is  froth  with  uncertain- 
ties and  difficulties.  The  latter  Is  partially  applicable  to  a seml- 
emplrlcal  approach. 

*6ood  scatterers  In  Infrared  and  also  absorbers. 
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To  illustrate,  consider  the  generation  of  a smoke  cloud  by  artil- 
lery shells.  In  Fig.  4(a)  we  show  the  progress  of  cloud  in  the  "x“  - 
direction"  of  a 105mm  WP  fired  statically  and  dynamically.  First  we 
see  that  the  dynamic  curve  is  significantly  different  than  for  the  static 
case  which  can  be  associated  with  the  dispersion  of  the  smoke  material. 

In  Fig.  4(b)  we  see  that  there  are  large  variations  in  the  near  surface 
dispersal  of  the  smoke  fill  material  and  smoke  cloud. 

Currently,  the  model  described  in  the  next  section  considers  the 
mortar  and  artillery  as  an  instantaneous  source  which  extrapolates  to  a 
diameter  of  9.1  meters  and  a height*  of  2.73  meters  at  time  zero.  The 
AMSAA/JTCG  model  defines  the  source  as  a function  of  fill  weight,  e.g., 
the  value  for  the  155mm  is  -9  meters  and  height  -2.7-5  meters  whereas 
for  lesser  weight  munitions  the  values  are  somewhat  less.  Although  the 
latter  numbers  are  consistent,  the  AMSAA  model  suggests  that  the  entire 
cloud,  assuming  a Gaussian  distribution  function.  Is  contained  within 
4x  the  latter  numerical  values  and  this  may  not  validate  the  data  avail- 
able from  static  tests  like  the  model  proposed  here.  Operationally  the 
live  shots  are  important.  This  fact  must  be  evaluated  In  the  model 
proposed  here. 

The  saving  feature  is  that  great  precision  Is  not  crucial;  only 
sound  estimates  for  operational  usage  are  necessary  to  establish  system 
performance  and  smoke  munition  effectiveness.  To  meet  these  requirements 
a smoke  model  will  be  described  in  the  following  sections. 

*For  non-white  phosphorus  munitions. 
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3.  MODEL  DEVELOPMENT 


3.1  GENERAL 

In  this  section  we  will  present  the  formulation  of  a smoke  model 
which  can  be  utilized  in  a combat  simulation  model.  The  intent  is  to 
address  all  features  including  any  shortcomings  since  the  model  described 
below  is  viewed  by  this  writer  as  a first  generation.  Considerable 
research  is  required  to  eliminate  uncertainties  in  the  following  items: 
input  data;  smoke  cloud  structure;  and  obscuring  aerosol  optical  proper- 
ties. 


3.2  SMOKE  CLOUD  STRUCTURE 

Since  we  are  concerned  with  the  time  history  of  a smoke  cloud 

produced  by  a military  munition(s)  we  have  drawn  upon  the  experimental 

results  reported  by  Rae  (10).  The  smoke  source  was  a single  generator  (oil 

fog)  which  is  a non-exothermic,  non-chemical  source.  For  wind  speeds 

greater  than  2m/sec  the  tests  showed  that  the  cloud  height  and  width 

can  be  defined  for  the  conditions  specified  in  Table  6 where  x is  the 
product  of  mean  time  wind  speed,  u (m/sec)  and  time,  t(sec). 

Table  6 

SMOKE  CLOUD  DIMENSIONS 


Turbulence  State 

Width  (m) 

■KOFJIXOMri 

High  (Lapse) 

9.14  + .419x 

2.73  + .137x 

Moderate 

(Slightly  unstable) 

9.14  + .238x 

2.73  + .073x 

Low  (Stable) 

9.14  + . 146x 

2.73  + .046x 

10Rae,  R.  W.,  "Smoke  Screenings  of  Vehicle  Routes, "(U)  CAORDE  Memo  19, 
June  1954  (Conf). 


Some  data  shows  that  the  slope  of  the  smoke  cloud  is  a constant  (z/x 
< 0.09)  for  wind  speeds  > 4.5  m/sec  and  approximately  0.17  at  -2  m/sec. 
The  results  in  Table  6 have  been  expanded  by  interpolation  to  cover  the 
range  of  Pasquill  stability  parameters  and,  are  listed  in  Table  7. 

Table  7 

CLOUD  SIZE  PARAMETERS  FOR  NON- EXOTHERMIC  MUNITIONS 


Email ra 

Width,  y (m) 

Height,  T (m)* 

Remarks 

Pasquill 

TTTiHJI 

A 

-3 

9.1  + .419x 

2.73  + .137x 

Sunny  day 

B 

-2 

9.1  + . 328x 

2.73  + .llx 

Day,  broken  clouds 

C 

-1 

9.1  + .238x 

2.73  + ,073x 

Overcast  day/night 

D 

0 

9.1  + .20x 

2.73  + ,06x 

Neutral 

E 

1 

9.1  + .18x 

2.73  + . 055x 

Evening/early  am 

F 

CM 

9.1  + .146x 

2.73  + .046x 

Evening/early  am 

*See  Appendix  G for  discussion  of  plume  height  for  Q <5  m/sec. 


It  should  be  noted  that  Soviet  measurements  with  a heated  smoke 
pot*  (AT=2°C) , for  a terrain  roughness, Z0  * 0.4  cm,  give  the  following 
cloud  slopes: 

tan  a * -.85  ( J-  ) + .07  ~ < 0 (unstable) 

L L (3.1) 

s .07  U - 0 (neutral-stable) 

where  a * tan"1  (z/X ) . This  result  may  be  comoared  to  the  values  in 
Table  7 with  the  aid  of  eqs.  (2.3)  and  (3.1). 

*“7  ” Tambient^ 
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When  Z0  * 0.4  cm  and  the  Pasquill  stability  class  is  A-B  the  value 
of  L_1  < - .07(m~^),  and  the  cloud  slope  = 0.11.  Since  the  ZQ  from 
Table  7 may  differ  somewhat,  the  value  of  0.11  is  a good  average  value 
(0.09  - 0.17)  for  category  B.  This  indicates  a good  correlation  between 
Pasquill-Turner  and  Monin-Obukov  stability  definitions. 

To  describe  the  geometric  shape  of  the  smoke  cloud  we  select  a 
semicone  for  the  near  ground,  point  source,  no n-exo thermic  smoke  munition. 
Thus,  the  volume  as  a function  of  time,  V(t),  is  given  by  the  expression 

V(t)  = f x(t)[-^-]2  (3.2) 

At  this  point  we  introduce  the  effect  of  a strong  exothermic  chemi- 
cal reacting  smoke  munition  like  white  or  red  phosphorus.  Here  the 
thermal  effect  overrides  the  momentum  exchange  mechanism  for  the  disper- 
sion of  the  smoke  aerosol  as  the  smoke  cloud  temperature  differential, 

AT,  is  non-negative  for  a long  period.  Further  if  (AT)  cloud  >0  the  cloud 
has  buoyancy  and  this  buoyancy  will  vary  with  the  atmospheric  gradient. 

Let  us  now  consider  the  thermal  history  of  the  smoke  cloud  and  estimate 
its  temporal  values. 

The  overall  heat  release  from  white  phosphorus  is  estimated  to  be 
800  kcal/kg.  Hence  the  Incremental  temperature  in  the  smoke  cloud, 
assuming  it  to  be  an  "air"  volume  is  given  by  the  simple  relation 

Qh  * mcpAT  (3.3) 

where  m,  Cp  and  qh  are  the  cloud  mass  (gms),  the  specific  heat  of 
air  and  the  quantity  of  heat  available  respectively 

For  a 10  meter  sphere  we  find,  where  m = p( density)  x V( volume), 

T - QH/cppV  - 10°C 
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At  later  times  AT  will  decrease  inversely  with  the  volume,  e.g.,  when 
the  radius  increases  to  ~20  meters  the  AT  ■*0.5°C,  and  therefore  some 
residual  buoyancy  remains  if  unstable  or  neutral  air  conditions  exist. 
On  the  other  hand  for  an  inversion  condition  (stable)  the  lift  effect 
would  stabilize.  We  will  consider  the  effect,  for  a more  practical 
case,  in  greater  detail. 

Consider  a continuous  source  then  the  heat  balance  equation  for 
combustion  is 


C * Qc  ■ S • pCpAT c 

(3.4) 

where 

C 

* concentration  (gm/m3) 

Qc 

3 heat  released  (cal  •m"^*sec"^) 

S 

3 gms  of  smoke  •m"^*sec"^ 

ATC 

* temperature  increase 

P 

3 smoke  density 

Similarly 

, the 

heat  released  by  reactive  processes 

including  atmospheric 

moisture 

reaction  is 

C • Qr  3 S • pCpAT r 

(3.5) 

where 

Qr 

3 T (pf  * rfNf  " pir3iNi)L 

(3.6) 

and. 

Pi 

3 particle  density 

ri 

3 particle  radius 

Ni 

3 no.  of  particles  formed/sec/meter 

L 

3 600  cal/gm 

Now  the  number  of  particles,  N remain  unchanged  and,  therefore 


I 

I 

E 

E 


e 

o 


0 


where 


Nf  = N.j  = v • S/(4npr3/3) 


v = fraction  of  hygroscopic  material 


(3.7) 


Hence  we  can  write  for  eq.  (3.6)  the  following 

Pfrf 

Qr  = SL  — r -1 
pfr1 


(3.8) 


The  effect  of  relative  humidity,  R%<100,  on  a hygroscopic  aerosol  can 
be  expressed  as 


TOO 


= 1 + 


3.2  x 10 
Tr< 


-5 


- 8.6 


m 


(3.9) 


where  m = initial  mass  in  grams 

W = molecular  wt  of  smoke  particle 

The  second  term  on  the  r.h.s.  is  evidently  small  and,  therefore  we 
have  for  r^, 

Q d m Amw  4 D 

(3.10) 


_ _ 8. 6m  _ n r 4tt  3,,,,,  R % 

rt-  77— =8.6  Tpr,/U(l-W) 


"0-  Toff) 

Substituting  eq.  (3.10)  into  (3.8)  gives 

Qr-'>  • s * lIt(7LV)-'] 

1 "TW 

and,  the  total  AT  is  given  in  eq.  (3.12), 


(3.11) 


AT 


PCr 


Qc 

-y-  n + v • L( 


(-  «!!■ ')] 


(3.12) 
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where 

Qc/S 

= 200 

n 

= combustion  efficiency,  smoke  wt/fill  weight 

CP 

= 0.24  cal/gm°C 

p 

= 1.29  x 103gm/m3 

V 

= ( Y-2 . 9 )/ Y 

Pf 

* 1.83  (phosphoric  acid) 

W 

= 98  (phosphoric  acid) 

ii 


If  we  substitute  the  appropriate  constants  then  for  R * 80%  and 
n = .9  we  obtain 

Li 

t‘7-  (.24K)?»«  l0»)  ('SO*  (.5)(600)  (^^-l)]*2.8CCC)  (3.13) 

For  an  uniform  concentration,  the  aT  is  uniform  throughout  the  smoke 
cloud  otherwise,  the  spatial  distribution  will  define  the  variation  in 
AT.  Since  AT  varies  with  relative  humidity  in  a non-linear  fashion,  we 
have  calculated  aT,  see  Fig.  5.  If  the  yield  curve  changes  upward  as 
suggested  by  Johnson  and  Forney,  the  aTs  will  be  differenct,  see  Fig.  5. 

Hence  if  C = 0.1  gm/m3  then  in  the  above  case  aT  = .28°C  a relative 
unimportant  value;  If  L * 100  and  CL  = 10  gm/m2  then  C = .1  gm/m3. 

To  a first  approximation  the  mean  concentration  can  be  determined 
by  dividing  the  smoke  weight  by  the  volume, 

(T  (x.y.z.t)  = Q/V  (3.14) 

i.e.,  we  assume  here  the  concentration  is  uniform  and  decreases  with 
time.  Although  this  approximation  suffers  from  at  least  two  shortcomings, 

0 
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Relative  Humdity  (*) 


Fig.  5.  Cloud  AT  as  a Function  of  Relative  Humidity 
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i.e.,  (1)  the  mean  concentration  has  a spatial  variation  and  (2)  local 
inhomogeneity  may  exist,  it  can  serve  to  provide  an  initial  value  for 
eq.  (3.12).  Further  the  shortcomings  will  be  examined  later. 

We  return  now  to  the  strong  exothermic  chemical  munition  and  deter 
mine  the  dimensions  of  the  resultant  smoke  cloud.  Here  the  thermal 
effect  overrides  the  momentum  exchange  mechanism  for  the  dispersion  of 
the  smoke  aerosols.  Since  the  cloud  temperature  differential  is  non- 
negative we  consider  the  plume  rise  to  be  defined  by  the  following 
expression 

Z = z + Ah  (3.15a) 

where  7 = momentum  rise  and. 

Ah  = rise  due  to  heat  addition 


For  Ah  we  accept  a modified  result  recommended  by  Briggs  (4)  when 


the  stability  category  is  unstable,  namely 

, 1/3  x 2/3 

Ah  = -LIE S (F1/6) 


(3.15b) 


where  F = 3.7  x 10  H and 
H = calories/sec* 


The  additional  factor  F ' accounts  for  the  difference  between  the  mean 
height  and  actual  cloud  top. 

When  a neutral  condition  exists  the  reconmended  expression  for 


Ah  is 


, ,Pl/2y*2/3r  y y 21[  y 1-2 

* | .4  +.64(ij+  2.2  (£j  j[l  +.8  (£*)j  (3.15c) 


*For  WP,  H = (800  cal/gm)  x WP  fill  wt(qm)  * 1 sec  burn  time. 


( 
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where  x*  = 0.52F'4hs'6 

hs  = source  height  = 1 meter 

The  maximum  rise  is  estimated  to  be  3x*s  x < 5x*. 

For  a stable  category  (3.15b)  is  good  to  x = 2.417/?^  where  ? is 
defined  by  eq.  I»  pll.  The  maximum  value  of  Ah  is  defined  by 

Ah(max)  = 2.9  (X)1/3  (3.15d) 

Since  z>y  the  smoke  cloud  shape  is  assumed  to  be  defined  as  a 
quarter  of  an  ellipsoid,  i.e.,  the  volume  is  given  by 

V(t)  - tt/6  • x *y  • Z (3.16) 

The  value  of  H for  various  munitions  is  uncertain;  that  is,  the 
rate  of  burn  and  conversion  are  unknown  and,  therefore  the  number  of 
calories  produced  per  sec  is  uncertain.  Theoretical  estimates  based 
upon  phosphorus  particle  size  and  rate  of  oxidization  would  suggest  =2 
sec  except  for  the  several  larger  fragments.  The  latter  implies  a 
dependence  on  explosive  charge  and  casing  thickness,  a fact  to  be 
established  by  experiment.*  Dolce  and  Metz  (11  ) state  "the  end  of  the 
weapon  (155  mm  shell)  phase  is  the  time  when  the  flowing  orange  hemis- 
phere becomes  a white  cloud,  which  is  generally  about  1 second  after 
burst."  The  conversion  of  the  initial  oxide  aerosol  to  phosphoric  acid 
droplets  occurs  at  distances  < 20  meters.  Hence  for  the  present,  until 

^ Dolce,  T.J.  and  Metz,  D.F.,  "An  Analysis  of  the  Smoke  Cloud  Data  from 
August  1975  Jefferson  Proving  Ground  Smoke  Test,"  AMSAA  Tech  Report 
201 , Sept  1977. 

*Private  communications  from  Dee,  Edgewood  Arsenal  (April  1978)  suggests 
for  a 81  mm  mortar  about  80'i  of  the  phosphorus  was  oxidized  in  < 2 sec 
and  the  remainder  < 10  secs.  Note:  the  source  was  ~5  meters  in  radius. 
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more  data  are  available  we  assume  the  conversion  is  =1  sec  for  all  WP 
munitions  irrespective  of  their  fill  weight  with  the  exception  of  the 
L8A1,  phosphorus  wicks  (WPW),  wedges,  etc.  The  latter  must  be  treated 
as  sources  with  a finite  duration  of  emission,  e.g.,  the  L8A1  burns  about 
30  secs  and  the  WPW  (2.75"  rocket  head)  burns  about  4 to  5 minutes. 

Further,  by  virtue  of  their  deployment  they  create  a linear  distributed 
source . 

Earlier  we  noted  that  the  temporal -spatial  distribution  defined 
by  a "Gaussian"  function  was  not  an  accurate  description  of  the  concen- 
tration in  a smoke  cloud  particularly  generated  by  military  smoke.  For 
the  vertical  distribution  we  are  aware  that  this  distribution  depends  on 
the  lapse  rate,  wind  speed  and  frictional  velocity,  u*.  The  horizontal 
distribution  of  the  concentration,  C(x),  depends  on  the  wind  speed,  surface 
roughness  and  a stability  length,  e.g.,  L,  the  Monin-Obukov  parameter. 

To  obtain  a qualitative  picture  of  the  (x,z)  distributions  assum- 
ing that  the  y distribution  is  uniform  we  examined  the  data  from  several 
experiments.  First  we  reduced  the  results  from  a Dugway  Test*  where  a 
linear  array  of  6-155mm  white  phosphorus  shells  were  statically  detonated.  The 
concentration  vs  x is  shown  in  Fig.  6 and,  the  extrapolation  to  x * 0 is 
given  in  Fig.  7. 

The  latter  corresponds  to  the  approximate  expression 

C(o,o,o,t)  = Coexp(-.0203t)  (3.17) 

where  C0  = 24.81  gm/m  . The  latter  is  consistent  with  initial  conditions. 

One  fact  that  this  indicates  is  that  the  mortar,  artillery,  etc. 
munitions  cannot  be  treated  as  a puff  since  the  residual  concentration 
is  “9  gm/m3  @ 60  secs.  Further,  the  effect  of  initial  dispersal  has  not 
been  considered. 

* Dugway  Trial,  DPI -002  T34,  Nov  19,  1977. 
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Fig.  7.  Temporal  Decay  of  Concentration  at  Source 
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The  concentration  C(x,o,o,t)  appears  to  be  defined  by  a relation 


K 

0 


C(x,o,o,t)  = C(o,o,o,t)  |l-Aef (3.18) 

where  A and  f(x-trt)  are  a numeric  and  a function  of  x respectively  and, 
C(x,o,o,t)-*-('-0.1 ) at  x = trt . This  dependence  is  clearly  related  to  a 
single  analyzed  event  where  the  stability  class  was  slightly  unstable. 
For  a similar  surface  situation  Kazanski  and  Monin  (12)  using  a contin- 
uous array  found  the  concentration  distribution  to  be  inversely  propor- 
tional to  (x/L)  which  will  agree  with  an  explosive  source  for  t»t0, 
i.e.,  quasi  - steady  state.  For  stable  conditions,  L increases  with 
stability  hence  the  decay  rate  will  increase  with  x.  Similarly,  for 
unstable  conditions  |1/L|  decreases  going  from  A-+0  stability  class  and 
therefore  the  decay  will  decrease  with  x as  we  go  from  stability  class 


A-0. 

l! 

For  the  vertical 

distribution  we  have  determined  for  a cloud  height 

n 

ij 

H that  the  dependence  can  be  categorized  as  follows: 

(unstable) 

C(z)  = C(x,o,o,t)  [1.25e‘Z/H-.25] 

0 

(neutral ) 

C(z)  = C(x,o,o,t)  [1.35e"Z/H-.35]  (3.19) 

0 

(stable) 

C(z)  = C(x,o,o,t)e1'52/H  for  Z < .2H 

(Z-.2H) 

* C(x,o,o,t)  [2.25e  H -.9]  forZs.2H 

0 

0 
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For  the  horizontal  distribution  we  suggest  the  following: 
a.  For  early  times  t s 20  secs 

C(x,o,o,t)  * C(o,o,o,o) [.8(1 .2  - x/ut)] 
where  x0  * ut 

— 

Kazanski,  A.B.,  and  Monin,  A.S.,  "The  Form  of  Smoke  Jets,"  Izv.  Akad. 
Nanka  SSSR  Seriya  Geoflzika,  1957,  p 1020. 
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b.  For  t>20 

C(x)  = C(o,  0,0,0)  (1  -x/xQ)  x/x0<.2 

C(x)  - C(o, 0,0,0) (1.1  -e(x/ut  _1))  x/x0  > .2 

For  a linear  array  equally  spaced  instantaneous  point  sources, 
e.g.,  n-munitions  spaced  "d"  meters  apart  where  d>y  (t=5)  they  can  be 
treated  as  individual  sources.  Fig.  8.  From  Fig.  8 we  note  the  overlap 

1 • 


Fig.  8.  Multi -Munition  Schematic 


occurs  at  time  t m (d  - 9.1 )/Atr  where  A is  the  constant  given  in  Table  6, 
y ■ 9.1  + Ax.  Beyond  this,  the  lateral  path  length  for  n = 2 will 
decrease  by  the  amount 

A*d«(x2  - x-|)/x-|  (3.20) 

where  x^  is  the  value  of  x at  time  of  intersection  and  x2  the  value  of  x 
at  the  point  of  interest.  For  n munitions  we  can  write  immediately 

A * (n-l)(d)(x2  -x1)/x1  (3.21) 

When  the  wind  speed  is  at  an  angle,  9,  appropriate  adjustment  is  required. 

The  overlap  zones  are  significant  for  radiation  transport  (optical 
properties)  laser  propagation  and  some  cases  of  backscatter.  However,  for 
C*L  computation  one  may  treat  the  sources  as  additives, £ C-jU. 
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For  a linear  continuous  source  of  length,  l,  where  Q grams  of  smoke 
material  of  yield  Y are  released,  the  quantity  of  smoke  produced  per  unit 
length  per  sec  is  (Q*Y)/(S.*t)  gms/m-sec. 

The  concentration  to  a first  approximation  is  defined  by  the  relation 

C(x,y,z,t)  = 2 (gm/m3)  (3.22) 

|u.2(x,t)(l+k) 


where  tr  = mean  speed  (m/sec) 

Z = smoke  cloud  height  @ "x"  (z  or  Z) 
k * incremental  length  of  smoke  width  @ "x" 

For  the  case  when  wind  speed  is  normal  to  the  smoke  line,  k * y/t,  and 
when  u is  at  angle,  0,  to  the  smoke  k * ysin0/t.  The  spatial  distribu- 
tions in  the  x and  z may  have  a form  similar  to  eqs.  (3.18)  and  (3.19). 

For  a distributed  source,  a combination  of  point  and/or  linear  arrays 
may  be  applied  with  appropriate  accounting  of  the  time  coordinates;  that 
is,  if  we  have  linear  arrays  at  x^Xj  and  xk  then  any  x is  (x-x^),  (x-xj) 
and  (x-xk)  and  the  relative  time  is  ti  * (x-x^)/ll,  etc.  A program  for 
distributed  steady  sources  like  smoke  pots  and  generators  is  given  in 
Appendix  0. 


When  a moving  ground  generator  is  utilized  a vectorial  representation 
is  applied;  that  is,  the  generator  velocity,  v,  is  added  vectorially  to  the 
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wind  velocity,  u,  hence  for  a headwind,  a line  source  is  produced  with  a 
resultant  mean  wind  speed  of  u + 7.  The  value  of  x is  found  by 


applying  the  expression,  x = [[u|  + |v|]t.  The  expansion  in  y and  z are 
defined  by  Table  7. 

When  the  generator  moves  at  an  angle  to  the  mean  wind  then  the 
resultant  speed  would  be  found  again  by  vector  addition. 


3.3  OPTICAL  CHARACTERISTICS 

The  overall  treatment  of  the  optical  characteristics  of  the  smoke 
cloud  and  its  interaction  with  optical /infrared  radiation  is  considered 
here.  Again  the  approach,  for  modelling  purposes,  is  to  provide  a simpli- 
fied version  of  radiation  transport  to  determine  the  smoke  effects  on 
optical  systems.  First  we  will  consider  the  aerosol,  then  the  transmit- 
tance and  finally  the  effects  on  systems. 

3.3.1  Aerosol  Properties 

The  notation  used  here  is  given  below 

a(X)  = spectral  cross  section  (m2)  or  (cm2) 

a(m)  = mass  coefficient  (m2/gm),  spectral  (X)  or  spectral 
band  (AX)  is  implied 

N = no.  of  particles/m3  or  per  cc 

a(X)  = No(X)  * volume  coefficient  (m~^) 

Subscripts 

e,s,a,  = extinction,  scattering,  absorption 
Definition 

extinction  = scattering  + absorption 

The  aerosols  produced  for  obscuring  purposes  vary  from  pure  scatterers 
(fog  oil)  to  particles  (phosphoric  acid)  which  are  highly  absorbing  In 
the  infrared.  Actually  the  obscuration  results  from  the  high  particle 
density  which  is  produced;  that  is,  for  a concentration  of  0.1  gm/m3  the 
number  density  is  of  the  order  of  105/cc  which  is  ~102to  103  greater  than 
that  found  in  clouds  or  haze.  This  fact  explains  the  larger  attenuation 
for  shorter  path  lengths  when  particle  sizes  are  comparable. 
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The  smoke  aerosol,  like  the  atmospheric  aerosol,  may  absorb  moisture 
(chlorides)  or  is  deliquescence  like  phosphoric  acid  and,  therefore  its 
optical  properties  will  depend  on  relative  humidity.  In  Fig.  9 the  mass 
extinction  coefficient  vs  wavelength  for  white  phosphorus  (WP)  is  given. 
From  Fig.  9 we  see  that  the  extinction  in  the  visible  is  significantly 
greater  than  in  the  infrared  and,  the  latter  varies  significantly  over 
the  region  of  several  microns  wide.  For  this  case,  integration  over  wave 
length  is  necessary;  that  is,  an  average  extinction  coefficient  cannot  be 
utilized.  In  regions  of  varying  absorption  it  has  been  established  that 
Beer's  law  is  not  valid.  Simply  speaking,  standard  methodology  should  be 
applied  namely,  to  integrate  the  product  of  (spectral  intensity  x system 
spectral  response  x spectral  transmittance)  over  the  wavelength  region  of 
interest.  In  the  event  a polydispersed  aerosol  mixture  is  produced  which 
is  the  usual  situation  then  the  mean  value  is  derived  by  integrating  the 
coefficient  as  a function  of  particle  size  over  the  size  distribution. 


Fig.  9.  Extinction  Coefficient  of  RP/WP 
;*ef  2,  Aoo  E.  E0-SP-7:002) 
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For  the  above  reasons  we  have  tended  here  to  depend  upon 
experimental  data  to  provide  the  mass  (nr/gm)  and  volume  (m-1)  optical 
parameters  (extinction,  etc)  for  the  different  smokes.  A precaution 
must  be  exercised  to  properly  account  for  relative  humidity:  the 
coefficients  will  increase  particularly  at  elevated  values. 

In  Appendix  E we  present  current  detailed  values  of  the  optical 
parameters  for  different  smokes.  Some  nominal  broadband  values  are  given 
in  Table  8. 

For  some  applications  there  is  a need  to  utilize  the  angular  dependence 
of  scattering.  For  small  particles  or,  when  the  particle  radius  < X/2ir 
(X  * wavelength  of  radiation)  then  the  scattering  Is  defined  by  (l+cos20), 
i.e.,  Rayleigh,  otherwise  the  scattering  is  highly  asymmetric  and  predomi- 
nantly in  the  forward  direction,  Mie  scattering.  To  obtain  the  scattering 
cross-section  we  integrate  over  all  solid  angles,  i.e., 

as  * /as(9)dfi, 

and. therefore  for  Rayleigh  scattering  the  scattered  intensity  flux  is 

equally  divided  between  the  forward  and  backward  direction.  For  Mie 

scattering  the  backscatter  is  relatively  small,  (102  - 103)  less  than  the 
forward  component.  By  convention,  one  normally  writes  os(e)  as  P(6)/4tt. 

In  Fig.  10  we  present  a result  for  phosphoric  acid  at  wavelengths 
3.39um  and  10.6um.  The  effect  of  wavelength  and  scattering  angle  are  self 
evident.  To  demonstrate  different  properties  we  also  plotted  the  value 
for  ZnCl2  (HC  smoke)  at  wavelength  of  10.5  um.* 

♦Private  communication  from  R.  Frickel,  CSL,  Aberdeen  Proving  Ground,  MO. 
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3.3.2  Radiation  Transport 

Here  we  will  present  the  necessary  formulations  to  determine  the 
effect  of  smoke  on  the  propagation  of  the  target  or  background  spectral 
intensity. 

3. 3. 2.1  Attenuation 

The  spectral  attenuation  is  defined  by  exp[-ae(A)R]  where  ae(A)  is 
the  spectral  extinction  coefficient  and  R Is  the  length  of  the  propagation 
path.  This  simple  expression  is  appropriate  when  (1)  scattering  is  omitted 
and  (2)  the  receiver  field-of-view  is  small.  For  a collimated  source  a 
first  order  correction  is  [1  + c^R  • $]  where  $ is  the  receiver  f-o-v/180°. 
Zuev  (13)  gives  a more  general  expression  for  single  scattering,  namely 
[1  + otgR  • D]  where  D is  a function  of  source  angle,  9,  receiver  angle,  <t>  and 
particle  size  parameter  p = 2irr/\,  see  Fig.  11. 


FIGURE  1 1.  Curv«  of  O'u  a function  of  receiver  aperture  :er  various  p 
and  half-cone  angles  9 of  the  source  emissions.  Solid  curves  pertain  to 
0 « a /S  and  dashed  curves  to 

For  smoke  the  mean  particle  size  is  -Z  urn  hence  at  10  urn,  p = 2irr/A  * 
1.256.  When  e » it/ 2 , a Lambertian  radiator,  and  <p  ^ 10  mrad  the  value  of  0 
is  insignificant.  Thus  for  a thermal  viewer  correction  for  single  scatter 
ing  is  unimportant.  If  A * 1.0  pm  and  a 10®  field  of  view  system  is  used 
then  the  correction  factor  is  sizeable. 

TS 

Zuev,  V.E.,  Propagation  of  Visible  and  Infrared  Radiation  in  the 
Atmosphere,  Halstead  Press,  John  Wiley  & Sons. 
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Clearly,  we  imply  that  the  ag(A)  is  treated  as  discussed  earlier. 
Hence  we  will  use  for  the  transmitted  radiation  the  following 

I = I0[l  +asR  •D(p,0,<O)]e‘aeR 


3. 3. 2. 2 Backscatter 

Here  we  must  define  two  concepts  of  backscatter.  First  the  back- 
scatter applicable  to  laser  radiation  corresponds  to  the  scattering  at 
180°,  P(tt).  This  can  be  obtained  from  plots  like  Fig.  10.  The  second 
concept  is  germane  to  solar  backscatter  or  other  sources,  e.g.,  flare 
illumination,  which  may  cause  glare  in  the  receiving  sensor,  see  Fig.  12. 

Under  the  assumption  that  single  scattering  is  valid,  i.e.,  when 
as/ae  s 0.8  and  multiple  scattering  can  be  neglected,  then  the  spectral 
steradiancy  (brightness)  is  defined  by 

T , cosi 


a- (A)  P(0) 

N\(e)  = Ns  (X )t-s  ( A ) • — j 7-v  ■ • t r— — 

A s'  ■ ' ap(A)  4it  cost  + cose 


(3.23) 


where  Ns  = source  steradiancy 

P(9)  * differential  scattering  cross  section/unit  area/steradian 

6 = tt  - (i  + e) , degrees 

i.e  * angle  of  incidence,  exitance 

T-j  * atmospheric  spectral  transmittance  from  source  to 
cloud  surface 

The  received  radiance  is  N\ie)x-re(A)  where  te(A)  is  the  spectral  trans- 
mittance from  the  cloud  surface  to  the  receiver.  For  an  uncoil imated 
source,  i and  9 will  include  a finite  set  of  values  and,  therefore 
N^(e)  will  require  integration  over  the  receiver  field-of-view.  For 
the  sun,  we  treat  it  as  a collimated  source. 


Receiver 


Laser 

Transmitter 


Receiver 


Fig.  12  Schematic  of  Backscatter  Geometries 


Again  under  the  assumption  of  single  scattering  the  albedo  of  a 
cloud  can  be  found  from  the  transport  equation  and  is  presented  in  Fig.  13; 
the  albedo  is  plotted  as  a function  of  optical  thickness,  aeR,  for  different 
values  of  as/ae.  For  a non-absorbing  medium,  as/ae  = 1,  the  albedo  approaches 
1 whereas  if  aa  is  finite  the  albedo  rapidly  approaches  an  asymptotic  value  <1 . 
The  latter  implies  that  multiple  scattering  effects  rapidly  decrease. 

With  this  brief  description  of  transmittance  and  backscatter,  we 
will  now  consider  contrast. 


3.3.3  Contrast 

The  contrast  concepts  discussed  here  cover  the  visible  and  infrared 
spectral  regions  where  the  possible  impact  of  thermal  effects  will  be 
examined.  Since  the  infrared  region  includes  thermal  viewers  where  the 
scene  information  is  displayed  such  systems  are  considered. 


3. 3. 3.1  Reflected  Radiation 

The  attenuation  of  target  and  background  radiation  is  only  necessary 
but  not  sufficient  to  describe  the  effect  on  an  target  acquisition  system 
of  an  obscuring  medium  in  a line  of  sight;  that  is,  the  observed  image 
contrast  is  a significant  parameter. 

By  definition  the  Inherent  contrast,  at  zero  range,  between  a target 
and  its  Immediate  background  is  expressed  here  by  the  equation,* 


Nt  - N 


B 


(3.24) 


♦This  definition  of  contrast  (Blackwell)  is  used  by  this  writer  as  part 
of  his  target  acquisition  routines  for  CARMOHETTE . 
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where 


Ny  = target  radiance,  and 
N0  = background  radiance 


The  observed  or  apparent  contrast  at  range,  R,  has  been  defined  as 

-1 


cR  • c0  n + 


1 -e'asR 


0.5(1  + Rg)e_cteR(RB) 


(3.25) 


where  Rg  * ground  reflectance  and 
Rg  * background  reflectance. 

A definition  for  apparent  contrast1  given  by  Middleton  (after  Duntley)(14)  is 


Cr*  co  l 1 + 4<e%R  -'> 


-1 


(3.26) 


where 


SK/Gd  * sky-ground  ratio. 


This  ratio  varies  with  ground  conditions,  cloud  cover,  and  scattering 
conditions  of  the  atmosphere.  Generally  for  a cloudy  sky,  the  ratio  is 
=8”1  and  0.2  8_1  for  a clear  sky  where  8 * Rg,  the  average  ground 
reflectance. 

For  daytime  conditions,  and  visual  systems,  CT,  the  threshold  contrast2  (50% 
probability  of  detection  level)  has  been  assumed  to  be  -2%  and,  the 
maximum  visual  range  (visibility)  V is  found  from  the  expression 


1 


A simple  derivation  can  be  given  as  follows;  if  G is  glare,  then 

<NT-NB).-“«R  6 a.R,-> 

R ' V*"  ° * 


2 

Sawyer,  K.F.,  Thorpe's  Dictionary  of  Applied  Chemistry,  4th  Ed.  (1950) 
reports  on  an  experiment  with  smoke.  The  objective  was  to  establish  the 
liminal  value  of  contrast.  Results  indicate  that  the  curve  had  the  same 
shape  as  Blackwell's  except  the  required  threshold  contrast  was  a factor  of 
two  greater. 

(14)  Middleton,  W.E.K.,  Vision  through  the  Atmosphere.  University  of 
Toronto  Press,  Canada  (1952). 
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C = e-ae*V  ; (3.27) 

T 

V 3 3.912/ae 

If  ae  is  defined  in  units  of  km~\  V is  in  kms. 

Duntley's  expression  omits  several  factors  such  as  sun  angle  and 
geometrical  relation  of  viewer  to  target.  Hruschke,  Rand  Corp.,  has 
developed  an  algorithm  for  the  determination  of  the  sky/ground  ratio  for 
different  conditions.  This  algorithm  is  now  presented. 

He  defines  the  ratio  by  an  expression 

SK 

g-p  F6(SG6  + SG£  ♦ SGv)/3  (3.28) 

where  SGg  3 exp(-l . 1 5 An  & -.75);  g * ground  albedo 

SG£  * exp(-3.4  sin  e + 2.7);  e 3 solar  elevation  (°) 

SGy  a exp(-0.5  An  V + 1.5);  V 3 visibility  (km) 

and,  <5  3 the  receiver  depression  angle. 

(1)  First,  for  cloud  cover  >4/8s,  e 3 25  . 

(2)  To  compute  F^  one  must  determine  a set  of  indices  for  different 
values  of  6,  e,  8,  and  V see  Table  9. 


Table  9.  - Values  of  Indices  jj,  j2»  J3.  j4  as  Functions  of 
5 » E » P » V 


5 (deg) 

j2  e (deg) 

j3  p 

J4  V (km) 

1 6 >82 

1 V>10 

2 50  >5  >30 

2 10  < e < 30 

2 . 5 > p > . 15 

2 10  > V 

3 30  > 5 > 12 

3 30  < e < 55 

3 . 15  > p 

4 12  > 5 

l I 1 
Li  1 

I 


(3)  With  the  established  values  of  the  indices  enter  Table  10  and 


find  F^. 

Table  10 

VALUES  OF  F6  AS  FUNCTION  OF  INDICES  j ^ j2>  j-j,  i* 


Jl 

^ 3 

^4 

F6 

*1 

j2 

j3 

^4 

F6 

jl 

^2 

^ 3 

^4 

F5 

Ji 

-*2 

•*3 

^ 4 

F6 

1 

1 

1 

1 

2 

2 

1 

1 

1 

1 

3 

1 

1 

1 

1.2 

4 

1 

1 

1 

1.2 

1 

1 

1 

2 

2 

2 

1 

1 

2 

1.2 

3 

1 

1 

2 

2 

4 

1 

1 

2 

0.8 

1 

1 

2 

1 

2 

2 

1 

2 

1 

1 

3 

1 

2 

l 

1.5 

4 

1 

2 

1 

1.5 

1 

1 

2 

2 

2 

2 

1 

2 

2 

1.2 

3 

1 

2 

2 

2.5 

4 

1 

2 

2 

1 

1 

1 

3 

1 

2 

2 

1 

3 

1 

1.2 

3 

1 

3 

1 

2.5 

4 

1 

3 

1 

2.5 

1 

1 

3 

2 

2 

2 

1 

3 

2 

2 

3 

1 

3 

2 

3.5 

4 

1 

3 

2 

1 

1 

2 

1 

1 

1.5 

2 

2 

1 

1 

1 

3 

2 

1 

1 

1.2 

4 

2 

1 

1 

1.2 

1 

2 

1 

2 

1.5 

2 

2 

1 

2 

1.2 

3 

2 

1 

2 

1.7 

4 

2 

1 

2 

0.8 

1 

2 

2 

1 

1.5 

2 

2 

2 

1 

1 

3 

2 

2 

1 

1.2 

4 

2 

2 

1 

1.2 

1 

2 

2 

2 

1.5 

2 

2 

2 

2 

1.2 

3 

2 

2 

2 

2 

4 

2 

2 

2 

0.8 

1 

2 

3 

1 

1.5 

2 

2 

3 

1 

1.2 

3 

2 

3 

1 

2 

4 

2 

3 

1 

2 

1 

2 

3 

2 

1.5 

2 

2 

3 

2 

1.5 

3 

2 

3 

2 

3 

4 

2 

3 

2 

1 

1 

3 

1 

1 

1.5 

2 

3 

1 

1 

1 

3 

3 

1 

1 

1 

4 

3 

1 

1 

1.2 

1 

3 

1 

2 

1.5 

2 

3 

1 

2 

1.2 

3 

3 

1 

2 

1.2 

4 

3 

1 

2 

0.7 

1 

3 

2 

1 

1.5 

2 

3 

2 

1 

1 

3 

3 

2 

1 

1.2 

4 

3 

2 

1 

1.2 

1 

3 

2 

2 

1.5 

2 

3 

2 

2 

1.2 

3 

3 

2 

2 

1.5 

4 

3 

2 

2 

0.8 

1 

3 

3 

1 

1.5 

2 

3 

3 

1 

1.2 

3 

3 

3 

1 

1.5 

4 

3 

3 

1 

1.5 

1 

3 

3 

2 

1.5 

2 

3 

3 

2 

1.5 

3 

3 

3 

2 

2 

4 

3 

3 

2 

0.8 

1 

4 

1 

1 

1 

2 

4 

1 

1 

1 

3 

4 

1 

1 

0.7 

4 

4 

1 

1 

1 

1 

4 

1 

2 

1 

2 

4 

1 

2 

1 

3 

4 

1 

2 

1 

4 

4 

1 

2 

0.5 

1 

4 

2 

1 

1 

2 

4 

2 

1 

1 

3 

4 

2 

1 

1 

4 

4 

2 

1 

1 

1 

4 

2 

2 

1 

2 

4 

2 

2 

1 

3 

4 

2 

2 

1.2 

4 

4 

2 

2 

0.7 

1 

4 

3 

1 

1 

2 

4 

3 

1 

1 

. 3 

4 

3 

1 

1 

4 

4 

3 

1 

1 

1 

4 

3 

2 

1 

2 

4 

3 

2 

1.2 

3 

4 

3 

•2 

1.5 

4 

4 

3 

2 

0.7 

aFor  82  > 5 > 50,  F6  = 1. 


At  this  point  we  consider  in  particular  the  smoke  layer  which  produces 
in  most  instances  "glare."  To  handle  this  induced  background  we  consider 
as  a first  approximation  that  (1)  the  ground  plane  reflectance  is  replaced 
by  the  smoke  layer  reflectance  and  (2)  the  Duntley  form  (eq.  3.26)  or  the 
Hruschke  algorithm  can  be  applied.  In  either  case,  the  albedo,  6,  of 
the  cloud  layer  is  defined  by  Fig.  13  and  when  the  visibility,  V,  is  required 
then  the  ambient  (unsmoked)  atmosphere  is  considered  as  the  propagation  medium. 
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Let  us  now  consider  the  contrast  modulation  which  is  given  by  the 
expression 

Nt  “ r 

CM  = Nt  + Nb  “ 2+C  (3-29* 

where  C,  the  contrast,  is  defined  above.  By  analogy  to  the  above  treatment 
in  "footnote,  p48",  we  can  write  eq.  (3.29)  as 

<CH)r  * «*><,['*  Jfjf  • e“eR  ] (3.30) 

Since  Q<C  and  C -1  is  a high  conspicuity  value  then  we  assume  (l/2+C)« 

0.4  and  G/Nr*^"1.  Hence  we  may  write  for  (C^)R 

(Cm)  R s (CM)o  [l+  (3.31) 

Let  B = .4  and  otR  * 1 then  the  contrast  modulation  is  reduced  to  about 
27%  of  the  original  value,  a serious  reduction  in  spatial  frequency 
content  making  the  recognition  task  difficult.  It  should  be  stated  that 
the  above  is  an  illustrative  example  and  estimate.  An  exact  calculation 
should  be  made  when  necessary. 

3.3.4  Thermal  Contrast 

Of  particular  interest  for  the  thermal  viewer  there  are  several 
factors.  In  addition  to  attenuation,  which  may  affect  picture  contrast. 
These  factors  are  the  following: 

1.  Channel  Crosstalk  - contrast  loss  due  to  scattering  within 
the  smoke  cloud  by  the  aerosol  constituents; 

2.  Temperature  - temporal  history  of  the  smoke  cloud  temperature; 

3.  Emission  - temporal  history  of  the  particulate  emission;  and 
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4.  Radiance  Modulation  - spatial  variations  of  radiance  due  to 
hot  spots  when  WP-like  material  is  used  and/or  concentration 
inhomogeneities  of  particulates. 

The  several  factors  are  germane  to  staring,  i.e.,  the  DC  type  thermal  viewers, 
whereas  the  factors  (1)  and  (4)  are  particularly  troublesome  to  the  current 
viewers  (AC  type).  We  will  now  consider  each  of  the  above  factors. 


1 . Channel  Crosstalk 

Let  <j>  be  the  system  instantaneous  field  of  view  (IFOV)  in 
milliradians  then  the  effect  of  smoke  could  be  such  that  a portion  of  the 
background  radiation  will  impinge  on  the  target  signal  resolution  element. 
Conversely,  some  target  radiation  can  fall  on  the  adjacent  cell,  see  Fig.  14. 


Background 


Target 


Fig.  14  Schematic  of  Scattering  Effect 

In  a similar  manner,  some  background  radiance  will  fall  on  the  target 
resolution  element.  We  describe  the  resultant  radiance  difference 
between  the  target  cell  (N^)  and  background  cell  (NB)  as 

AN  * [NS (0 ) + NB(o1)]  - [ NB ( 0 ) + NS(^)]  (3.32) 

n 
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where  N-j (4>-j ) is  the  scattering  contribution.  Rewriting  eq.  (3.32)  we 
obtain 

( Np  (<J>-i  ) - Nq(4>'>  ) j 

AN  = (NS(0)  - NB(0)]]l  - j 

If  we  identify  (AT)p  = AN(0),  the  effect  of  scattering  is 


We  will  now  tacitly  assume  that  single  scattering  will  occur  and 
this  is  reasonable.  That  is,  for  a WP  screen  (.02  gm/mJ)  the  attenuation 
length  is  -0.144  kms  and,  for  a diesel  fog  dispersal  (.005  gm/m3)  the 
attenuation  length  is  -4.34  kms.  The  latter  two  numbers  bracket  the 
overall  spectrum  of  values.  Since  the  likelihood  is  small  that  a trans- 
lucent smoke  layer  will  be  thicker  than  *150  meters  the  assumption  is 
reasonable. 

We  can  now  calculate  Ns(<j>-|),  Ng(d>i ) by  the  use  of  simple  diffraction 
theory.  The  reason  Is  that  <j>i  *0.50  milliradian. 

The  scattered  flux,  F(e),  is  given  by  the  equation 
F(0)  * F(0)e_aeL  • 0(e) 
where  for  a monodispersed  aerosol 


the  parameters  are 

N » particle  concentration/km3 
p * ird/X  (d  * diameter) 

z * pe  (kilometer-radian) 

J-j  a Bessel  Function  of  the  first  kind 
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].  Substituting  the  appropriate  values  for  an  instantaneous  field  of  view 

= 0.25  mrad  we  obtain  the  values  of  AN(4>)/AN(0)  =0.03  or,  the  reduced  AT 

~0.97ATo.  Hence  crosstalk  can  be  dismissed  as  a problem. 

n 

2.  Temperature 

In  Section  3. 2we  presented  an  exposition  of  the  thermal  history 
of  the  smoke  cloud  and,  showed  AT  <2.8°C  x concentration  for  a relative 
humidity  <80%.  These  values  of  AT  are  not  detrimental  with  the  possible 
exception  of  a high  performance  DC  system  with  MRT  -.01. 

Lj 

3.  Particle  Emission 

The  conversion  of  chemical  compositions  to  smoke  is  normally 
an  exothermic  reaction  producing  particulates  that  may  range  in  the  200°C 
to  500°C.  In  the  specific  case  of  phosphorus  there  is  an  exothermic 

reaction  that  converts  the  oxide  particle  to  phosphoric  acid  droplets.  The 

(J 

initial  conversion  energy  is  transferred  to  the  early  phase  of  the  smoke 
cloud  whereas  for  phosphoric  acid  the  droplet  remains  hot  -several  hundred 
degrees  for  a finite  time.  The  question  is  how  long  will  it  be  hot. 

The  temperature  history  at  r.  along  radius  of  a spherical  particle 
(radius  = a)  at  an  uniform  temperature,  T0,  is  given  by  (15)  as 

Tn  2KaTn  * 

T (r’t}  3 KTT  * -IT9-  l*  g n x f(an»r) 

where  an = roots  of  tan  aa  = 3aa/3+ka2a2 

K = ratio  of  medium  heat  capacity  to  particle  heat  capacity  and, 
k * thermal  conductivity 

When  kt/a2  >0.16  then  T(r)  * T0/K  + l.  Substituting  appropriate  values 
15 

Carslaw,  H.S.,  and  Jaeger,  J.C.,  Conduction  of  Heat  in  Solid,  Oxford,  1947. 
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into  the  former  we  obtain  t = 1.6*10”6  sec  which  is  due  to  the  small  particle 
radius,  a,  being  only  -1  pm.  To  estimate  K we  consider  a smoke  concentration 
~1  gm/m3  then  K-300  and  the  particles  will  have  equilibrated  with  the  air 
in  the  smoke  cloud  since  we  have  shown  earlier  that  AT  =2.8°C  at  a concen- 
tration of  1 gm/m3.  This  appears  to  be  consistent  with  observations;  that  is, 
except  for  early  times  (C-L  large)  the  particle  emission  is  not  significant. 
Hence  for  C*L>  system  threshold,  particle  emission  is  a non-contributory 
factor  to  system  performance  degradation.* 

4.  Radiance  Modulation 

For  the  present  we  will  only  discuss  the  mechanism  in  a quali- 
tative manner  as  a detail  analysis  is  beyond  the  scope  of  this  study.  In 
the  above  we  have  considered  the  smoke  cloud  density  to  be  uniform  or 
locally  isotropic  except  in  the  regions  of  overlap.  Similarly  in  the 
gaussian  approach  no  local  homogeneities  exist.  Hence  neither  approach 
provides  for  local  small  scale  eddies  arising  from  atmospheric  shear  or 
terrain  features.  In  the  latter  case  the  effects  are  localized  and  are 
comparable  to  or~2Xthe  size  of  the  terrain  feature.  If  the  thermal  viewer 
system  resolution  provides  0.3  meter  at  the  smoke  cloud  then  the  inhomo- 
geneities will  be  confined  to  the  ground  level.  Other  inhomogeneities  at 
practical  times  would  be  comparable  to  a major  fraction  of  the  cloud 
dimension  or,  significantly  larger  than  thermal  resolution  elements. 

Hence,  local  hot  spots  from  munition  fragments,  burning  wicks 
or  spatially  distributed  smoke  pots  may  cause  serious  problems.  Also 
whenever  an  overlap  condition  occurs  there  may  be  occasions  that  these 

*N0T£:  This  is  not  to  be  confused  with  (1)  hot  "chunks"  at  t*0  secs  or 
(2)  hot  smoking  chunks  which  can  cause  target  blockage. 
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overlaps  may  give  rise  to  AT(s)  comparable  to  a target's  AT  when  the  target 
is  detected  through  the  smoke  cloud.  In  this  general  category  the  vertical 
stratification  may  cause  a problem  in  the  vicinity  of  the  peak  (for  stable  or 
inversion  conditions).  Similarly,  the  observed  fluctuations  in  lateral  concen- 
tration due  to  axial  wind  fluctuations  may  under  special  circumstances  lead 
to  a signal  modulation  problem  in  time  and  space  coordinate. 

To  summarize,  we  suggest  that  operational  modelled  experiments  be 
performed  to  examine  conclusively  if  circumstances  exist  where  the  smoke 
cloud  contains  regions  that  can  appear  as  clutter.  For  modelling  we  propose 
to  use  the  several  expressions  for  contrast  and  transmission  germane  to  the 
proposed  model . 

First,  the  initial  obscuration  is  induced  by  the  cloud  opacity, 
namely  high  attenuation.  The  transmission  over  the  range,  R,  between 
target  and  observer  is  given  by 

T ■ (1  + a|-4D)  exp  -[(a*ir  (R-4)  + a*  £)]  (3.34) 

where  J l and  a*  are  the  smoke  path  length  and  linear  extinction  coefficient 
of  smoke  respectively.  The  value  of  l is  a function  of  observer-target 
orientation  with  respect  to  plume  axis  and,  is  discussed  in  Appendix  F. 


The  second  effect  is  contrast  reduction  of  D.C.  systems  by  the 
glare  or  backscatter  from  the  plume. 


Hence  we  propose  as  a first  order  approximation  to  utilize  the 
following  expression  for  the  perceived  contrast  through  smoke,  Crs. 

Crs  s Cr  -Ts  • | l+a/6(easR*  - 1)  -1  (3.35) 
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where 


Cr  3 contrast  without  smoke 


T$  = transmittance  through  smoke 

a = constant  (1,  0.2  for  non-sunny  and  sunny  conditions 
respectively) 

R*  = range  from  cloud  edge 

a$  * smoke  cloud  scattering  coefficient 

To  obtain  a precise  backscatter  estimate  the  phase  function  may 

be  used. 

In  Volume  2 we  have  assembled  the  model  details  into  a formal 
structure. 


4.  VALIDATION 


In  this  section  we  present  the  results  of  a validation  of  the  material 
described  in  the  earlier  section;  that  is,  on  the  basis  of  the  agreement 
between  the  model  formulated  earlier  and  results  from  several  field  tests 
we  conclude  that  the  proposed  model  can  provide  reasonable  prediction  of 
the  smoke  cloud  history  and  the  obscuring  effects  of  smoke.  We  now 
present  the  results  on  cloud  rise,  cloud  width  and  then  actual  obscuration 
effects. 

4.1  Cloud  Rise 

A significant  parameter  of  cloud  dynamics  is  the  rate  of  cloud  rise 
as  a function  of  munition  fill  material,  namely  exothermic  (phosphorus) 
and  other  smokes.  Specifically,  we  evaluate  the  impact  of  rate  of  heat 
release  on  cloud  height,  since  the  cloud  height  affects  cloud  volume 
magnitude. 

For  the  purpose  of  validation  we  selected  the  results  of  a number  of 
munitions  tested  at  the  Jefferson  Proving  Ground  (10).  The  comparative 
analysis  is  given  below. 

(a)  60mm  Mortar 

Fill  - .75  # white  phosphorus,  burn  time  = 1 sec 
Relative  humidity  - 62% 

Pasquill  Cat  C 

Mean  Windspeed  - 2 m/sec 

With  the  above  conditions,  we  first  applied  the  appropriate  expression 
for  the  cloud  rise,  z(t)  as  given  in  Table  7.  This  result  is  shown  in 
Fig.  15.  We  then  calculated  the  increment.  Ah,  due  to  the  use  of  WP, 
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using  eq.  (3.15).  The  total  cloud  height,  Z(t),  is  compared  with  the 
observed  values  where  the  effect  of  the  variation  in  IT  is  included. 

The  agreement  between  the  computed  and  observed  results  is  excellent. 

This  and  subsequent  results  support  the  need  and  method  of  including 
thermal  effects  into  the  formalism  proposed  here. 

(b)  81mm  Mortar 

Fill  - 1.75  # white  phosphorus,  burn  time  = 1 sec 
Relative  humidity  - 92 % 

Pasquill  Category  = A 
Mean  Windspeed  = .5-1.0  m/sec 

In  Fig.  16,  we  present  the  comparative  result  for  the  81mm  mortar.  The 
agreement  is  satisfactory  at  later  times.  Uncertainty  in  the  wind  speed 
makes  a precise  correlation  difficult. 


Fig.  16  Plume  Rise  vs  Time:  81mm 
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(c)  105mm 

Fill  - 3.83  # white  phosphorus,  burn  time  = 1 sec 
Relative  humidity  * 92% 

Pasquill  category  = A/B 
Mean  windspeed  - 1.0-2  m/sec 


The  results  are  given  in  Fig.  17  and  show  good  agreement. 
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(d)  155mm 

Fill  = 25.84  # 155  HC 
Relative  Humidity:  92% 

Pasqulll  Category:  B/C 
Mean  winds peed:  1.2-2. 5 m/sec 

In  Fig.  18  we  present  the  comparative  results  and  again  with  good 
agreement. 

An  item  of  significance  is  the  fact  that  the  cloud  approaches  an 
asymptotic  value  as  discussed  in  eqs.  (3.15). 

Hence  we  can  state  that  the  plume  rise  algorithms  including  the 
effect  of  cloud  buoyancy  is  valid. 


o — obs 
■ — • calc. 
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Fig.  18  Plume  Rise  vs  Time:  155mm  HC 


4.2  CLOUD  WIDTH  AND  OBSCURATION  EFFECTS 

(1)  The  effects  of  the  smoke  cloud  on  electrooptical  systems  are,  as 
noted  earlier,  dependent  upon  the  product  (concentration  x smoke  path 
length).  Hence  the  ability  to  calculate  the  concentration  along  the  line- 
of-sight  is  basic  to  any  model.  In  the  validation  tests  to  follow,  the 
line-of-slght  is  established  a'priori  since  it  is  the  line  which  contains 
the  particle  density  sampling  devices. 

(2)  The  cloud  width  is  defined  in  Table  7 for  single  smoke  munitions. 

To  determine  the  total  cloud  width.L.as  a function  time  for  several  colllnear 
munitions  one  may  add  to  the  initial  width  at  time  zero  the  width,  y( t), 
from  Table  7,  see  Fig.  19. 

jy/2- 

l 

} a/2. 

Fig.  19.  Schematic  for  Cloud  Width  Determination 


(3)  Here  three  tests  for  validation  are  considered;  the  first  two  at 
White  Sands  Missile  Range  and  last  at  Dugway  Proving  Ground.  One  should 
note  that  the  test  ranges  are  relatively  low  humidity  (semi-arid)  environ- 
ments and,  the  ground  is  relatively  smooth.  For  example  at  Dugway  the 
relative  humidity  ~35X.  This  fact  tends  to  suppress  the  thermal  effects. 
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The  first  test  was  to  compare  the  results  obtained  from  the  static 
detonation  of  5 - 4.2"  WP  rounds.  The  input  data  were  relative  humidity  - 
31%,wind  speed  = 3.4  m/sec,  wind  direction  250°  and  Pasquill  Category  B-C. 


The  geometrical  configuration  and  surface  cloud  history  are  shown  in 
Fig.  20.  The  calculated  width  Is  ~80m  whereas  the  observed  width  is  ~83  m. 

In  Fig.  21  we  present  a comparison  of  the  calculated  and  measured 
(CxL)  as  a function  of  time  where  the  rive  rounds  are  treated  as  five 
individual  events  neglecting  the  effect  of  overlap  zones  and  considered 
the  volume  to  be  ellipsoidal.  The  agreement  is  fair  since  we  have  omitted 
fluctuations  in  the  wind  speed,  however  the  calculated  mean  (CxL)  would 
have  accurately  predicted  the  effect  on  a thermal  system.  In  addition, 
an  estimate  of  the  atmospheric  transmission  t in  the  8-12pm  region  cor- 
responding to  CxL  values  Is  depicted  In  Fig.  21. 

Let  us  now  consider  the  detail  shape  of  the  CxL  curve.  From  later 
results  we  find  that  the  time  mean  wind  speed,  u-,  provides  the  mean  con- 
centration; however,  any  fluctuations  may  be  significant.  Since  we  are 
interested  in  fluctuations  on  the  order  of  seconds  with  sampling  times 
~1  sec  and  these  data  are  unavailable  we  can  only  suggest  that  in  part 
these  fluctuations  are  responsible.  The  reason  is  that  the  overlapping 
may  contribute  additional  turbulence. 

In  Fig.  22  we  have  reproduced  the  observed  temporal  history  of  con- 
centration at  the  three  aerosol  samplers.  We  note  that  the  two  exterior 
samplers  are  similar  whereas  the  center  sampler's  behavior  is  entirely 
different  but  consistent  with  the  overlap  zones.  The  observed  gaps  in 
concentration  may  be  due  to  wind  fluctuations. 
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Fig.  21  C*L  History  for  5-4.2”  WP  Mortars 


One  should  note  that  the  above  estimate  is  of  practical  value; 
however,  an  improved  methodology  is  required  to  treat  the  overlap  region. 
The  reason  is  that  the  outer  layers  are  optical  thinner  than  the  inner 
core  and,  the  effects  on  radiation  transfer  may  be  different  than  for  an 
uniform  concentration  along  a 1 ine-of-sight,  e.g.,  when  considering 
backscatter. 

Validation  (#2)  - WSMR  Test  12,  WPW-2,  21  July  1977 

The  second  validation  is  more  complex  and  represents  a non-steady 
state  continuous  line  array.  Here  13-  2.75"  rocket  heads  were  detonated 
in  a static  mode.  Each  head  contains  12  white  phosphorus  wicks  (WPW) 
and  the  total  weight  of  white  phosphorus/head  is  1.682  kg.  The  meteoro- 
logical conditions  were  relative  humidity  - 36%  and  mean  wind  speed  - 3.7 
m/sec.  The  Pasquill  stability  category  was  stated  to  be  B. 

In  Fig.  23  is  the  geometric  layout  of  the  smoke  material.  We  assume 
in  the  absence  of  any  distribution  data,  that  the  source  is  three  lines 
where  the  centerline  (~75  m from  sampler  line)  contains  50%  of  source 
strength  and  50%  is  equally  distributed  + 20  m from  the  centerline. 

First  we  note  that  the  effect  of  the  WP  is  insignificant  since  each 
wick  has  only  '-140  gms  of  slow  burning  phosphorus.  Hence  the  width 
increase  and  height  of  the  cloud  are  defined  by  Table  7.  The  values  at 
15,  20,  and  25  secs,  arrival  times  from  the  three  lines  to  sampling  line, 
are  y * 27.1,  33.4  and  39.5  and  z * 8.8,  10.98  and  13.2  m respectively. 


Smoke 


Sampling 

Line 


Fig.  23  Sensor/Smoke  Layout 


From  the  above  the  maximum  width  is  l + y * (120  + 39.4)m  -160m.  The 
observed  maximum  width  -150  meters. 

We  now  compute  the  concentration  in  steps,  i.e.,  from  each  source 
line  separately.  The  concentration  at  the  sampling  line  is  given  by 


C(x,t) 


1 

£u  Z(x,t)(l+k) 


(3.36) 


where  Q 
Y 
t 
l 
U 

z(x,t) 

k 


* smoke  material  wt  (gms) 
a yield 

* burn  time  (sec) 

* length  of  smoke  material  (m) 
a mean  speed  (m/sec) 

* smoke  cloud  height  at  x 

* [smoke  width  (y)  at  x]/: 
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Substituting  the  values  for  constants  in  (3.36)  @ 14.9  secs  to  obtain 


C(x,t) 


[5.46  x 1 0 3 ) ( 3 ) 
[240) (120) 


-JT  (3.7)(8.8)(1 .22) 


= .055  g/m3  @ 14.9  secs. 

At  20.2  secs,  the  smoke  from  the  2nd  line  arrives  and  the  smoke  from  the 
3rd  line  arrives  at  25.6  secs.  The  individual  concentrations  are  found 
directly  from  (3.36) 

C(x,20)  = .055  x 2 x = .088  g/m3. 


C(x,25.6)  = .055  x 1 x ^ = .037  g/m3 

The  parameter  (C*L)  = concentration  x smoke  thickness  and  the  individual 
contributions  are 

(1)  C( x, 1 5 ) = .055  x 147  = 8.10  gm/m2 

(2)  C(x,20)  = .088  x 153  = 13.5  g/m2  or,  @20  sec  = 21. 6 g/m2 

(3)  C(x,25)  = .037  x 160  = 5.9  g/m2  or  @25  sec = 27.5  g/m2 

At  later  times  we  assert  that  the  mean  concentration  drops  as  M(t)/V(t) 
and  therefore,  we  obtain  the  following  values  at  40  and  120  secs,  namely 
~17.54  g/m2  and  -9.4  gm/m2,  see  Fig.  24. 

The  oscillatory  character  of  the  "C*L"  curve  clearly  depends  on 
the  turbulence  spectra  of  atmosphere  during  the  test;  that  is  u"  is  the 
mean  speed  whereas  the  actual  speed  is  u = H + u'  where  u'  is  the  fluctua- 
tions. To  appropriately  find  the  variance  of  the  wind" speed  we  need 
sampling  times  ~1  sec  for  a sampling  period  of  ~10  mins  i.e.  t0  + 5 mins. 
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An  examination  of  the  wind  speed  at  2 and  7 meters  sampled  every 
60  secs  indicated  that  a minimum  occurred  at  +150  secs  and  a peak  occurred 
at  '•210  secs  an  excursion  of  -1.75  to  2 in  C.  The  standard  deviation  is 
-30$;  however,  the  excursion  of  a factor  of  2 can  reduce  or  increase  the 
C*L  by  a factor  of  2.  The  reason  is  that  the  concentration  «tr"^  . The 
mean  value  and  the  standard  deviation  are  shown  in  Fig.  24. 

Further  validation  has  been  obtained  from  a comparison  between 
predicted  and  observed  system  outage  but  omitted  for  classification 
reasons. 


Fig.  24  C*L  History  for  13  - 2.75"  Heads: 

(-  ) measured;  (-•-)  calculated 
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Validation  (#3)  - Dugway  Test  DPI-002-T-34,  19  Nov  1977 

Here  we  consider  the  static  test  firing  of  six  (6)  155  mm  (M110E2) 
white  phosphorus  shells  spaced  3.7  meters  apart.  The  input  data  are  as 
follows:  mean  speed  - 4.9  m/sec  at  20°  to  horizontal  axis;  relative 
humidity  - 35%;  and  Pasquill  category  - D. 

Before  presenting  any  calculational  results  let  us  briefly  discuss 
the  experimental  details.  The  close  proximity  of  such  large  shells  raises 
the  question  on  the  interactive  effects  and  the  dispersal  of  the  smoke 
source  material.  Further  the  size  of  the  phosphorus  particles  is  uncer- 
tain and,  in  turn  the  burn  time. 

The  test  data  include  the  cloud  size  at  time  of  detonation:  height 
=7  m;  width  =25  m;  and  length  =22  m.  Hence,  the  phosphorus  material  was 
dispersed  over  the  area  25  x 22  m2.  Interestingly,  since  the  original 
width  is  18.5  m the  additional  linear  spread  was  only  ~6.5  m.  If  we 
now  consider  an  individual  round  and  apply  the  appropriate  Ah  equation* 
we  find  x*  = 4.4  m and  Ah  -*•  23  m in  several  seconds.  Thereafter  the  plume 
will  rise  gently  by  momentum  exchange. 

Since  the  phosphorus  cannot  burn  long  enough  to  be  treated  as  a 
continuous  source  we  assume  that  source  is  approximately  instantaneous. 

The  (C*L)  result  is  shown  in  Fig.  25  for  the  Sampling  Line  at 
x ~76  m.  For  comparison  we  show  the  measured  C*L  at  x -1 37  m and  we 
note  that  the  two  C*L  values  as  a function  of  time  are  approximately 
similar.  The  latter  is  an  inherent  assumption  of  the  model,  namely 
uniform  concentration. 

*For  category  D,  eq.  (3.15c ) . 
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Fi9-  25  C*L  History  of  6-155mm  WP  Shells 
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Appendix  A 

MUNITIONS  CHARACTERISTICS 

I 


Mb /(Obsolete)  FS  4. 59 

Hi/,  Hb/Al(SIU.C)  UP  4.06 
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Table  A2 

Some  Smoke  Generators  and  Pots 


i. 

11 

Ii 


Smoke 

Type 

Fill 

Fill  Wt/Flow  Rate 

Smoke  Generator 

M3 

Fog  Oil 

25-50  gals/hr 

Smoke  Pot 

AN-M7 

Fog  Oil  (SGF) 

— 

AN-M7A1 

Fog  Oil  (SGF) 

— 

Smoke  Pot 

Ml 

HC 

10.0 

ABC-M5 

HC 

3.0 

Air  Smoke  Pot 

M4A2 

Table  A3 
U.K.  Grenades 


Smoke 

Screen 

Type 

1 

WHSm 

■jPQj 

Ural 

Source 

Smoke 

— 

■ral 

Fill 

L5A1 ,2,3,4 

70  (grd) 

55-75 

l l 

White 

8 

L7A1 

70  (grd) 

30-50 

Point 

ground 

Green 

XL7E2 

120  (grd) 

30-50 

Point 

ground 

Green 

XL8E1 

(Curtain) 

L8A1 

30 

60-90 

8 m air, 
ground 

2 

RP+ 

butyl  rubber 

XL18E1 
(for  L8A1 
practice) 

120 

10  m air, 
ground 

w 

silica 

XL19E1 

no 

8-10 

Green 
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DISCUSSION  OF  SMOKE  YIELD 
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DISCUSSIOM  OF  SMOKE  YIELD 


The  yield  of  smoke/ obscurants , particular  whi.te  phosphorus,  is  uncer- 
tain. Johnson  and  Forney ^ published  a set  of  data  which  has  been  reduced 
here  in  algebraic  form: 

Y(WP)  = 3.8  + .003  (RH%  - 10)1,67 
Y(ZnCl2)  = 1 + .051  (RH%  - 5)'85 
Y(A1C13/H2S04)  = 2.8  + .016  (RH%  - 20)1 *25 
Y(Fog  Oil)  = 1 


L! 


Other  data  indicate  that  yield  for  white  phosphorus  can  be  defined  by 
the  expression,  eq.  2.2. 

Y (WP)  = 2.9  + .003  (RH%  - 40)1'67 

The  other  obscurant  yields  have  adjusted  accordingly.  This  difference 
is  significant  and  can  seriously  alter  the  validation  of  a model;  that 
is,  the  spatio-temporal  history  of  the  aerosol  concentration  and  radia- 
tion transfer. 

Another  impact  of  the  yield  is  the  effect  on  AT  (eq  3.13).  Since  the 
parameter  v relates  the  liquid  fraction  to  the  total  mass  we  may  also 
write  v = (Y-k)/Y.  If  k = 2.9  or  3.8  the  value  of  v will  change  accord- 
ingly and  similarly  aT,  see  Fig.  5. 


^Johnson,  M.C.  and  Forney,  P.D.,  "The  Effectiveness  of  Obscuring  Smokes," 
ORG,  Edgewood  Arsenal,  1972  (unpublished). 
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Appendix  C 


A STABILITY  CLASSIFICATION  BASED  ON 
HOURLY  AIRPORT  OBSERVATIONS 


A STABILITY  CLASSIFICATION  BASED  ON 
HOURLY  AIRPORT  OBSERVATIONS* 

This  system  of  classifying  stability  on  an  hourly  basis  for  research 
in  air  pollution  is  based  upon  work  accomplished  by  Dr.  F. 

Pasquill  of  the  British  Meteorological  Office  (see  reference  1 of 
text) . Stability  near  the  ground  is  dependent  primarily  upon 
net  radiation  and  wind  speed.  Without  the  Influence  of  clouds, 
insolation  (incoming  radiation)  during  the  day  is  dependent 
upon  solar  altitude;  which  is  a function  of  time  of  day  and 
time  of  year.  When  clouds  exist,  their  cover  and  thickness 
decrease  Incoming  and  outgoing  radiation.  In  this  system 
insolation  is  estimated  by  solar  altitude  and  modified  for 
existing  conditions  of  total  cloud  cover  and  cloud  celling 
height.  At  night,  estimates  of  outgoing  radiation  are  made 
by  considering  cloud  cover.  This  stability  classification 
system  has  been  made  completely  objective  so  that  an  electronic 
computer  can  be  used  to  compute  stability  classes.  The  stability 
classes  are  as  follows:  (A)  Extremely  Unstable,  (B)  Unstable, 

(C)  Slightly  Unstable,  (D)  Neutral,  (E)  Slightly  Stable, 

(F)  Stable,  (G)  Extremely  Stable.  Table  A-l  gives  the 
stability  class  as  a function  of  wind  speed  and  net  radiation. 

The  net  radiation  index  ranges  from  4,  highest  positive  net 
radiation  (directed  toward  ttje  ground),  to  -2,  highest  negative 
net  radiation  (directed  away  from  the  earth) . Instability 
occurs  with  high  positive  net  radiation  and  low  wind  speed, 
stability  with  high  negative  net  radiation  and  light  winds, 
and  neutral  conditions  with  cloudy  skies  or  high  wind  speeds . 

The  net  radiation  index  used  with  wind  speed  to  obtain  stability 
class  is  determined  by  the  following  procedure: 

* The  following  explanation  of  the  Pasquill  Stability  classification 
has  been  extracted  from  an  article  by  D.  Bruce  Turner  in  the 
February  1964  Journal  of  Applied  Meteorology 
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(1)  If  the  total  cloud  cover  is  10/10  and  the  ceiling  is  less 
than  7000  feet,  use  net  radiation  index  equal  to  0 (whether  day 

or  night). 

(2)  For  nighttime  (night  is  defined  as  the  period  from  one 
hour  before  sunset  to  one  hour  after  sunrise) : 

a.  If  total  cloud  cover  <4/10,  use  net  radiation  index 

equal  to  -2. 

b.  If  total  cloud  cover  >4/10,  use  net  radiation  index 

equal  to  -1. 

(3)  For  daytime: 

a.  Determine  the  Insolation  class  number  as  a function 
of  solar  altitude  from  Table  A-2. 

b.  If  total  cloud  cover  <5/10,  use  the  net  radiation  index 
in  Table  A-l  corresponding  to  the  insolation  class  number. 

c.  If  cloud  cover  >5/10,  modify  the  insolation  class  number 
by  following  these  six  steps: 

(1)  Ceiling  <7000  ft.,  subtract  2. 

(2)  Ceiling  £7000  ft.  but  < 16000  ft.,  subtract  1. 

(3)  Total  cloud  cover  equal  10/10,  subtract  1.  (This 
will  only  apply  to  ceilings  £7000  ft.  since  cases 
with  10/10  coverage  below  7000  ft.  are  considered 
in  item  1 above) . 

(4)  If  insolation  class  number  has  not  been  modified 
by  steps  (1),  (2),  or  (3)  above,  assume  modified 
class  number  equal  to  insolation  class  number. 

(5)  If  modified  Insolation  class  number  is  less  than 
1,  let  it  equal  1. 

(6)  Use  the  net  radiation  index  in  Table  A-l 
corresponding  to  the  modified  Insolation  class 
number. 


I 
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TABLE  A-l  STABILITY  CLASS  AS  A FUNCTION  OF  NET  RADIATION  AND  WIND  SPEED 


WIND  SPEED  NET  RADIATION  INDEX 

(KNOTS)  4 3 2 1 0 -1  -2 


0,  1 
2,  3 
A,  5 
6 
7 

8,  9 
10 
11 
> 12 


TABLE  A- 2 INSOLATION  AS  A FUNCTION  OF  SOLAR  ALTITUDE 


SOLAR  ALTITUDE  INSOLATION 

(a)  INSOLATION  CLASS  NUMBER 


60*  < a 

Strong 

4 

35*  < a 

60 

< Moderate 

3 

15*  < a 

35 

S Slight 

2 

a i 15° 

Weak 

1 
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DESCRIPTION  OF  ATDL  COMPUTER  MODEL  FOR 
DISPERSION  FROM  MULTIPLE  SOURCES  » 


Steven  R.  Hanna 

Air  Resources  Atmospheric  Turbulence  and  Diffusion  Laboratory 
National  Oceanic  and  Atmospheric  Administration 
Oak  Ridge,  Tennessee 


In  order  to  estimate  surface  concentrations  of  pollutants  in  regions  con* 
taining  many  sources,  it  is  convenient  to  place  the  emissions  into  two 
categories:  area  sources  (e.g.,  home  heating  or  small  industries)  and  point 
sources  (ej..  power  plants  or  large  industries).  In  previous  reports,1  ’ 2 simple 
methods  of  calculating  concentrations  due  to  area  sources  were  derived  and 
verified.  We  have  devised  a straightforward  computer  model  for  accounting 
for  point  sources  as  well  as  area  sources.  This  model  has  been  applied  to  many 
regions,  including  Rockwood-Harriman,  Tennessee,2  and  Knoxville,  Tennes- 
see.4 It  is  the  purpose  of  this  paper  to  outline  the  ATDL  computer  program 
in  detail  so  that  others  may  use  it. 


CALCULATION  OF  SURFACE  CONCENTRATIONS  DUE  TO  AREA 
SOURCES 

Most  source  emission  inventories  distribute  the  area  sources  in  an  urban 
region  into  a square  grid  pattern,  with  grid'distances  of  about  S km,  where  it 
is  assumed  that  the  area  source  strength  in  any  grid  square  is  uniform  across 
tfiat  square.  Gifford’s2  “reciprocal  plume”  concept  is  employed  in  order  to 
estimate  the  surface  concentration  yfp% /mJ)  due  taarea  sources  upwind  of 
the  receptor  point: 


where  QA0rg/m2sec)  is  source  strength,  D(m)  is  the  distance  to  the  edge  of 
the  urban  area,  U(m/sec)  is  the  wind  speed,  and  oz(m)  is  the  vertical  disper- 
sion parameter.  Smith4  suggests  the  following  empirical  form  for  az: 

a,  * ax'2,  (x  in  meters)  (2) 


*Noll,  K.E.,  (ed) , Industrial  Air  Pollution  Control,  Ann  Arbor  Science 
Publishers  Inc.,  1975. 


Sunny  day 

a ■ .40 

b ■ SI 

Cloudy  day 

a - .15 

b m .75  (also  applies  to  yearly  averages) 

Night 

a ■ .06 

b * .71 

If  the  receptor  is  at  the  center  of  grid  block  “0,”  with  grid  distance  Ax, 
and  the  wind  blows  in  only  one  direction,  then  equation  1 can  be  written  as 
the  summation  over  grid  squares  upwind  of  the  receptor  square: 


X 


(Ax/2)‘-b 
Ua(l -b) 


♦XQ* 


i*l 


The  source  strengths  Qaq,  Qa«  ’ ^A2 *PP*y  to  the  grid  square  in  which 

the  receptor  is  located,  the  grid  square  upwind  of  the  receptor  square,  and  so 
on,  respectively.  The  integration  is  arbitrarily  terminated  after  four  grid 
squares.  If  the  frequency  with  which  the  wind  blows  from  the  16  major 
directions  is  known,  then  equation  3 becomes  the  double  summation: 

, r ±4  t4 

. . fTl  (W  b I 0.  (0.0)  ♦ I I Q*(i.i)  f (id) 

X fT  UaO'ST  I A 

(4) 

[(2r+l)l'b  - (2r  - l)l  b 


where  r is  the  number  of  grid  blocks  that  square  (i.j)  is  away  from  the  central 
receptor  square.  For  example,  r equals  one  for  the  “ring"  of  eight  squares 
adjacent  to  the  central  receptor  square.  When  the  16  point  wind  direction 
frequency  distribution  is  input  to  the  program,  the  parameters  f (ij)  and 


(2r+l)I*b  • (2r-l)1*b  are  calculated  within  the  program,  using  techniques 


explained  by  Gifford  and  Hanna.1 

In  the  program,  therefore,  each  element  C(lj)  of  a 9 by  9 matrix  of 
coefficients  is  multiplied  by  the  corresponding  element  QA(lj)  of  the  source 
matrix  to  obtain  the  concentration  x in  the  central  square.  Thus,  it  is  neces- 
sary to  expand  the  source  matrix  by  4 squares  on  all  sides,  entering  small, 
background  source  strengths  in  these  squares,  so  that  the  9 by  9 coefficient 
matrix  will  “fit”  over  the  source  grid  along  the  edges  of  the  given  urban  area. 

These  techniques  are  straightforward  and  do  not  necessarily  require  a 
digital  computer.  We  have  often  made  them  using  a slide  rule  or  desk  calcula- 


-PA* 
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tor.  However,  the  point  source  calculations,  described  in  the  next  section,  are 
sufficiently  complicated  that  a digital  computer  should  be  used. 

CALCULATION  OF  SURFACE  CONCENTRATIONS  DUE  TO  POINT 
SOURCES 

The  basic  Gaussian  plume  formula  described  by,  for  example,  Slade7  is 
used  to  estimate  surface  concentrations  due  to  point  sources.  Fora  16  point 
wind  direction  frequency  distribution,  the  surface  concentration  x of  a pollu- 
tant emitted  at  strength  Qp  Gig/sec)  is  given  by  the  formula: 


o,rU2l 
1 16 


where  H(m)  is  the  effective  source  height,  r(m)  is  the  distance  of  the  receptor 
point  from  the  source,  and  f is  the  frequency  with  which  the  wind  blows 
towards  the  sector  of  interest.  (This  formula  is  equivalent  to  Slade’s7  equa- 
tion 3.144.)  The  effective  source  height  H is  the  sum  of  the  stack  height  h, 
and  the  plume  rise  h_,  which  is  calculated  using  the  formula  proposed  by 
Briggs:*  K 

hp  - 2.9|w0R02  [T^  -Teo]  / [U(*J-  + .01*  C/m)] | »/3  (6> 


where  g,  T,  wQ,  and  R are  the  acceleration  of  gravity,  absolute  temperature, 
initial  plume  vertical  speed,  and  stack  radius,  respectively.  The  subscripts  jj 
and  e refer  to  plume  and  environment  variables,  respectively. 

Equation  5 can  be  used  in  order  to  estimate  the  average  concentration  in 
the  grid  square  in  which  the  point  source  is  located.  The  source  is  assumed  to 
be  located  in  the  center  of  the  square,  and  the  average  concentration  in  the 
circle  with  radius  Ax/2  is  assumed  to  equal  the  average  in  the  square  with  side 
Ax.  Also  since  there  is  only  one  wind  direction  sector  for  the  source  square, 
the  fraction  2rr/16  in  equation  5 becomes  2n,and  the  frequency  f equals  one. 


source  square 


2a2r2b 
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This  equation  is  integrated  numerically  by  the  computer,  given  the  input 
parameters  Ax,  Q , a,  b,  H,  and  U.  It  is  found  that  the  avenge  concentration 
due  to  a ground  level  (H  ■ 0)  point  source  located  in  the  middle  of  the  grid 
square  is  equal  to  the  avenge  concentration  due  to  a uniform  area  source  in 
that  grid  square,  provided  that  Qp  equals  QA  (Ax)2 ; Le.,  the  total  emissions 
in  the  square  are  equal.  The  average  concentration  x in  grid  squares  other 
than  the  source  grid  square  is  assumed  to  equal  the  concentration  in  the 
center  of  the  square  calculated  using  equation  5.  Gifford's*  reciprocal  plume 
concept  is  then  used  to  estimate  the  total  surface  concentration  in  a grid 
square  due  to  contributions  from  the  given  distribution  of  point  sources  by 
superimposing  a 9 by  9 matrix  over  the  source  grid,  and  multiplying  term  by 
term.  This  is  done  in  the  program  by  classes  of  effective  stack  height,  H. 

PROGRAM  DESCRIPTION 

Input 

The  computer  program,  written  in  FORTRAN  IV  for  use  on  an  IBM 
360/6S  computer,  performs  the  above  calculations  using  the  following  input 
parameters.  Maximum  permissible  dimensions  and  units  are  given. 

Card  1;  (2F1 0.2,415) 

DX:  grid  distance  (m) 

BX:  rural  or  background  source  strength  Org/m2  sec) 

NR:  number  of  rows,  including  an  extra  4 rows  on  the  top  A bottom 

NC:  number  of  columns,  including  an  extra  4 columns  on  either  side 

Ng:  number  of  corrections  to  area  source  strengths  (if  no  corrections, 
put  Np  » 1 . If  no  area  sources,  put  Np  * 0.) 

NH:  number  of  effective  source  height  classes 
Card  2;  (1615);  NN:  number  of  stability  classes 

Card  3;  (8F10J);  RX(  10):  parameters  a in  az  • axb 
Card  4;  (8F10.2);  PX(10):  parameters  b in  at  ■ axb 
Card  S;  (8FI0.2);  U:  wind  speed  (m/sec) 

Card  6;  (3F10.2);  F(16):  wind  direction  frequency  distribution  begin- 

ning with  NNE  and  going  clockwise 
Card  7;  (20A4);  mTl : description  of  area  sources 

Card  8;  (20A4);  ITIT2:  description  of  point  sources 

Card  9;  (20A4);  OTTO:  description  of  concentration  patterns 

Card  10;  (1615);  ID(30):  row  number  of  area  source  correction.  If  no 
correction,  ID(1)  * 1 

Card  11;  (1615);  JD(30):  column  number  of  area  source  correction.  If 

no  correction,  JD(l)  ■ 1 

Carl  12;  (8F10.2);  ST(30):  area  source  corrections  Org/m1  sec).  If  no 

correction,  ST(1)  ■ .0 

Card  13;  (8F10.2);  S(30,30):  area  source  strengths  Org/m2sec) 
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Card  14;  (8F10.2);  H(20):  effective  source  heights  (m) 

Card  IS;  (1615);  NP(20):  number  of  point  sources  in  each  source  height 
class 

Card  16;  | f (16IS);1C(10»30>:  row  number  of  point  source 

Card  17;  > < (16IS);JC(  10,30):  column  number  of  point  source 

Card  18;  ) ((8FI0.2);QEC(10,30):  point  source  strength  (pg/sec) 

If  no  point  sources,  put  NH  ■ 1,H(1)«  1,NP(1)*  1,  IC(1,1)* 
1,JC(1.1)-1,QEC(1  ,l)-0.0 


Output 

Maximum  permissible  dimensions  are  given. 

1.  (8F10.2);  F(16):  wind  direction  frequency  distribution,  beginning  with 

NNE  and  going  clockwise 

2.  (8F10.2):  FU(9,9):  9 by  9 matrix  of  direction  frequency  distribution, 

divided  by  wind  speed 

3.  NR  (number  of  rows),  NC  (number  of  columns),  N$  (number  of  area 

source  corrections),  NH  (number  of  effective  source  height  classes) 

4.  DX  (grid  distance  [m]),  BX (rural  source  strength  [pg/m2secj),  U (wind 

speed  [m/sec]) 

5.  R,  B:  current  values  in  o2  * Rx®,  applying  to  the  following  numbers  6-13 

6.  C(9,9):  matrix  of  coefficients  for  area  sources 

7.  A(30,30):  concentrations  (pg/m3)  due  to  area  sources 

8.  H(20):  current  effective  source  height,  applying  to  numben  9-1 1 

9.  IC(10,30),  JC(  10,30),  QEC( 1 0,30):  current  location  and  magnitude  of 

point  sources 

10.  SS(9,9):  matrix  of  coefficients  for  current  point  sources 

1 1 . AA(30,30):  concentrations  (pg/m3)  due  to  current  point  sources 

12.  AC(30,30):  concentrations  (pg/m3)  due  to  all  sources 

13.  RAT(30,30):  ratio  of  concentration  due  to  area  sources  to  concentre* 

tion  due  to  all  sources 

EXAMPLE  OF  OUTPUT  OF  PROGRAM 

Consider  the  town  of  Moonshine,  Tennessee  (population  603),  a square 
town  which  has  been  divided  into  a square  grid  system  consisting  of  four  S 
km  by  S km  grids.  Each  grid  square  has  an  area  source  emission  of  1 
pg/m1  sec.  A point  source,  with  effective  source  height,  H * 30m,  and  source 
strength  10T  pg/sec,  is  in  the  center  of  the  southwest  block.  The  wind  blows 
with  equal  frequency  from  ail  directions; Le.,  F(l)  ■ F(2) ...  * F(16)  *0.67. 
Assume  ihat  the  wind  speed  U equals  S m/sec.  Consider  only  the  stability 
class  applicable  to  average  yearly  conditions:  NN  * 1,  RX(1)  * .15, 
PX(I)  * .75. 
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Other  input  parameters  are  then:  DX  • 5000.,  BX  * .0,  NR  = 10,  NC  » 10, 
Np  - 1 , NH  - 1 , 1D(1)  - 1 . JD(l)  - 1 . ST(1)  - .0,  S(5  fi)  - S(5 .6)  « S(6,5)  - 
S(6,6)  - 1..  NP  - 1,  IC(1 ,1)  - 6,  JC(1 ,1)  - 5.  QEC(1 ,1)  ■ 107 . 

The  predicted  surface  concentrations  due  to  area  and  point  sources  are 
listed  in  Table  3.  These  numbers  can  be  used  by  anyone  to  check  whether  or 
not  they  are  using  the  program  correctly. 


Ifcble  3:  Source  Data  and  Predicted  Surface  Concentrations  at  Moonshine, 


Surface 

Concentration 


Surface 

Concentration 


Area 

Point 

Due  to 

Due  to 

Grid 

Source 

Source 

Area  Sources 

Area  Sources 

Square 

pg/sec  m1 

107  Mg/aec 

/ig/m’ 

jig/m* 

NW 

1.0 

0.0 

34.64 

0.58 

NE 

1.0 

0.0 

34.64 

0.33 

SW 

1.0 

1.0 

34.64 

3.70 

SE 

1.0 

0.0 

34.64 

0.58 

APPENDIX:  PROGRAM  LISTING 

Btasastoa  stio.JO)  . ratio. ]0)  .ctio.ioi  .aoo.ioi  .Qttjo.ioi  .ctio.ioi , 
iss ijo.io) .ntio.ioi  .ictio.ioi.xctic.ioi .jctio.io.arcno. joi . roti 
20)  .JO  (30)  .ST  (100)  .IX(10).IS(10|«  ZTXT1  (20)  .ITIT2(20)  .ITITJ  (20)  .1 (20 
i) . ur (20)  .us(io.io)  .aiT(io.iO)  .rut) 

100  rOIRIT  (2F10.  2.  *XS) 

ios  ro»ai?(erio.2) 
no  roast?  (20i«) 
in  rotaiTtittS) 

allots. ioO)  Ds.ai.aa.ic.io.ia 
do  loo  x-i.aa 
oo  ioo  j-i.ic 
ioo  six.ji-ix 
mots. i*ii  ae 
mots. ios)  did). x«i.ii) 
snots. ios)  trxti). x-i.n> 
mots. ios)  o 
mots,  ios)  inn  .x-i.i*) 

ItlOIS. 110)  XTXT1 
snots. 110)  ITXT2 
■nets.1101  x Tin 

ItlDlS.  1*1)  (Xt(X).X-I.IO) 

1110  IS.  1*1)  (JO  (J|  « J-1 . 10) 
snots.  1051  (ST  |J|  .J-1. 10) 
jr-ec-s 
xj-ea-s 

mots. ics)  ns(x.j)  .j-s.jx)  .x-s.xi) 
mots,  ios)  ii tx)  .x-1. an 

*110 IS.  1*1)  tltlX)  .X-1.ll) 

DO  1*0  X-1. IS 
air-iitx) 

anots.i*})  txctx.j) .j-i.xii) 
mots. i»ii  uctx.j)  .j-i.sii) 

1X10(5.1051  IQICIT.JI  .J-1. Ill) 
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miS  FAQB  IS  BEST  QUALITY  PRACTICAL 
mOM  OOPX  fTJRKISHED  TO  DDC 


MO  COITIIOI 

OtX«Ilt.704|  (fill.  1-1. 141 
704  PO lilt (21.000X00  OOXX/(Ori0.4|| 

00  700  1*1.0 
00  700  3*1.0 
7oo  roix.3i*.o 
1-1. /» 
ro  is. si  •( 

00  700  1*1.4 

rois-x.s«xi*ri2i*i 
ro  is»z.s*xi*nti*o 
ro  is*z.s-xi*riioi  *c 
ro  is-x.5-ii*rn*i  *t 
rois-x.si*rnti*o 
ro  is.s»zi*ri4i  *o 
rois»x.si*rtoi  *t 

700  rois.s-xt*ni2i*o 

ro  i4.si«roi4.si*i*(riit  ♦rnsn 
ro  is.oi*ro(S.(i»i*iri)i*risi  i 
ro  io.si«roit.si  ♦«*iri7i  *r ioi  i 
rois,4i*rois.4i*i*(riiii*rMJii 
DO  701  Z*1. 2 
ron*z.oi*riii*t 
roi4.o*zi*riJi*i 
ro  it.t«xi  *nsi  *t 
roit*x.ti-ri7i  •« 
ro  io*;.4)*rioi*t 
ro  it.i»zi*riiii*i 
roi4.i«zi*riu»*t 

701  roii*z.4t*rrisi*c 
ro  n.7i*r  m *r 

ro  u.oi*ri3i*i 
ro  i7.oi*risi  *r 
ro  io.7i  *r  i7i*« 
ro  10.31  *r«oi*i 
ro  i7.ii-rnii*i 
roi3.ii*r<i3i*i 
roii.3i*rnsi*t 

MZTXIt.O)  liroiZ.JI  .J*1.0l  .Z*1.0I 

0 roioit  ixx.nozotot  mo  sotto  in  mooooct  n?axi/iorio.sn 

(IZTtlt.4701  OI.OC.IO.OO 

470  rOOOAT 1 21.0000.  IC0S*.ZS.3t. 12000.  COLOOOS*.  ZS.  31.3000*.  ZS.  30,  2100 

io  trrtCT  soooct  otx-.xsi 

OtZTtlt.lOl  01.01.0 

io  roooiT i2z. oooi  zo  a *.ri. o.3i. 2IOOOOU  soosci  STieooroM.rs.  1.3Z.1 
noozoo  sroio*.no.3.soo/sKi 
00  100  Z*1.0I 

00  too  j*1.0C 
too  Ql  IZ.JI  *.  0 
00  If 7 1*1.00 
Z*  ZD IOI 
3*30101 

107  SIZ.3I*STI0I 

OIZTIIt.2001  ZTZT1.  IISIZ.3I  .>  1 .OCI  . Z-1 . 00| 

200  roriuTi2i.20i*/iiorio.2)i 

DO  102  ZZ'1,00 
0*00  IZZI 
t*OX  IZZ) 

oizTXit.noi  o.o 

no  roioitizx. *300(0  sot  or  roott  no  rnioxTitt  roo  tzeoi  >/3x.2to*.r 
IS.  3.3X.2O0*.  PS.  31 
zriooi  ito.ito.iso 
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The  appended  curves  are  reproduced  from  the  following  two 

reports : 

(1)  Milham,  M.E.  et  al , "New  Findings  on  the  Nature 
of  WP/RP  Smokes,"  USARRADCOM  TRARSCL-TR-77067, 
July  1977:  Figs.  D-l , D-2. 


(2)  Milham,  M.,  "A  Catalog  of  Optical  Extinction 
Data  for  Various  Aerosols/Smokes,"  Edgewood 
Arsenal  Special  Pub.,  ED-SP-77002;  June  1976. 


Experimental  Low  Relative  Humidity  WP/RP  Spectrum 
(J-  to  Sfim  region) 


Expcrimmt.il  Low  Relative  Humidity  tvp.'RP  Spectrum 
♦ 7-  to  1 4-ttm  region  i 
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Appendix  F 
SMOKE  PATH  LENGTH 


For  illustrative  purposes,  we  present  here,  a brief  discussion 
of  the  geometry  of  the  path  through  a smoke  cloud. 

Consistent  with  the  model  of  a right  cone  we  will  derive  some 
expressions  for  the  path  length,  L,  through  the  smoke  plume. 

We  consider  the  smoke  source  to  be  the  origin  of  coordinates 
(0,0,0)  and  consider  the  observer  and  target  coordinates  to  be 
(x0.yo»o)  and  (xt*yt’°)  respectively.  See  Fig.  E-l 


/(vyt.o)  -t 


Fig.  E-l  Geometrical  Arrangement 


From  Fig.  E-l  we  can  write  the  following  for  the  coordinates  of 


the  segment  through  the  plume 

X1  ■ (y0  + y0 tan#)/ (Cot  # + tan  0) 
y-j  3 -tan  0 • 

x2  3 (y©  + xQtan#)/(Cot  # + tan  0) 
y2  3 tan  0 3 (9.1  + Ax)/2x 

and,  its  length,  L,  is 

L * 2(yQ  + x0  Cot  <p) (Cot  <(>  +0)/(Cot20  • Cot2#  - 1 )** 


where  tan  0 = (9.1  + Ax)/2x 

It  should  be  noted  that  in  the  limit,  lim  Cot  0 3 2/A. 

x-*» 

When  OP-T  is  normal  to  the  x-axis  then  L 3 2xQ  Cot  0. 

The  above  would  also  apply  to  the  WP  smoke  cloud  geometry,  i.e., 
it  is  the  near  ground  case.  More  generally  the  path  length  is 

L * 2 - [t»tT >-h(0P ) J 2 

where  Ah  » |Z  (observer)  - Z (target) | 

R 3 slant  range  to  the  target 
rg(t)  « half  width  of  the  plume  (y/2) 


f 


NOTE:  If  several  munitions  are  detonated  nearby  simultaneously  in  the 
same  zone,  i.e.,  overlap  on  expansion,  then  the  plume  width  and 
concentration  can  be  treated  as  additive  provided  the  spatio 
temporal  history  is  accounted  for  correctly. 
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Appendix  G 
PLUME  HEIGHT 


In  Table  7 we  have  defined  the  vertical  rise  as  2.93  + A*x  where 
A depends  upon  Pasquill  category.  For  modelling  purposes  we  have  proposed 
a fixed  value  for  each  category.  Actually  the  parameter  "A"  varies 
continuously  with  mean  wind  speed,  u,  for  values  <5  m/sec.  Although 
smoke  normally  is  not  discharged  in  u<2  m/sec  because  of  the  horizontal 
coverage  is  low,  for  completeness  we  present  the  value  of  A for  0.5 <u 
<5  m/sec.  The  latter  will  cover  stability  categories  A,  B,  and  C. 
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